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Abstract 

A classification of plienomenologically interesting supersymmetric extensions of the Standard- 
Model with a U{1)r symmetry is presented. Some of these are consistent with subsets of leptonic 
or baryonic "R-parity violating" (RPV) operators, thereby providing a natural motivation for them. 
We then focus on a particular class of models in which the U{1)r symmetry coincides with lepton 
number when restricted to the SM sector. In this case, the extension of lepton number to the su- 
perpartners is "non-standard" , implying, in particular, the existence of the leptonic RPV operators 
LLE^ and LQD'^, and a vacuum structure where one of the left-handed sneutrinos acquires a signif- 
icant vacuum-expectation-value, while not being constrained by neutrino mass bounds. The model 
can be naturally consistent with bounds from electroweak precision measurements and flavor-changing 
processes. It can also easily accommodate the recently measured Higgs mass due to the existence of 
a scalar triplet that couples to the Higgs with an order one coupling, with only moderate fine-tuning. 
The phenomenology is rather rich and distinctive, with features such as heavy-but-natural Dirac gaug- 
inos, relaxed bounds on squarks, resonant slepton/sneutrino production, lepto-quark signals, as well as 
an interesting connection to neutrino physics arising from i?-breaking. The broad qualitative features 
are discussed in this paper, with a more detailed phenomenological study carried out in a companion 
paper [l]. 
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1 Introduction 



We are extremely fortunate to live in a data-rich era of particle physics. The discovery of a Higgs- 
like particle [2||3] is indeed a monumental achievement of the LHC. Measuring the properties of this 
particle in detail is now one of the most important experimental tasks. On the other hand, the mass 
of this particle (~ 125 GeV), as well as null results for beyond the Standard-Model (SM) physics so 
far, have started challenging our simple expectations for physics beyond the SM|^ In particular, the 
Minimal Supersymmetric Standard Model (MSSM), which is the leading candidate for beyond-the- 
Standard Model (BSM) physics, is being significantly constrained. For example, the bounds on colored 
superpartners are well over a TeV in the "bulk" of parameter space, where pair production of squarks 
and/or gluinos is followed by decay into quarks or gluons and the lightest supersymmetric particle 
(LSP), giving rise to a jets plus^T signature. 

How can one interpret the Higgs-like discovery at ~125 GeV, and the null results for new physics 
so far? Our discussion will be within the supersymmetric paradigm for concreteness. The current data 
seems to suggest two rather different approaches to beyond-the-SM (BSM) physics. One viewpoint is 
that a Higgs near 125 GeV and the absence of superpartners so far, can be explained in the MSSM 
with heavy (~ 10 — 100 TeV) scalars. These models are therefore electroweak-tuned. Models of this 
type can arise in some versions of [4||5]. In fact, there exist top-down frameworks which predict heavy 
scalars in the 10-100 TeV range |6||7] and, therefore, a Higgs mass near the observed value js]; the big 
hierarchy problem of dynamically generating the (~10) TeV scale from the Planck scale is solved in 
such a framework, but a little hierarchy remains. 

However, the other viewpoint is that it is premature to completely abandon electroweak naturalnes^ 
at this stage, for it is still possible to imagine models in which the bounds on superpartners are evaded 
in a natural manner, allowing electroweak-natural (at least to a large degree) models. At the same 
time, various mechanisms exist which can give rise to a Higgs-like particle at around 125 GeV without 
introducing excessive tuning. With this philosophy in mind, some approaches to (SUSY) BSM physics 
have generated renewed interest, such as: i) Models in which the first and second generation squarks 
are rather heavy, but those of the third generation are light for some reason, as in |9[]15|, ii) Models 



which have a rather special spectrum, such as a compressed one 16 18 , or a stealth one [19 ,20 , iii) 



Models which do not give rise to signatures with appreciable missing energy, the prototypical example 
being that of baryonic i?-parity violating (RPV) models (see 21 for a review), and iv) Models in 
which the production cross-section is small even for light first and second generation squarks, as in 



examples with a relatively heavy Dirac gluino 22 , 23 



In 24 , a model was proposed that shares some features of the RPV models and models with Dirac 
gluinos, mentioned above. The defining feature of this model was the existence of a t/(l)i? symmetry 
that was identified with one of the lepton numbers and the role of the sneutrino as the down-type Higgs. 
In 25 the R-symmetry was generalized to a global lepton number to allow for viable neutrino masses 
and mixings. Ref. |[26] studied the case where the U{1)r symmetry is identified with the baryon number. 
In this paper we perform a classification of phenomenologically interesting i?-symmetric models, and 



-"^For concreteness, we interpret the ~ 125 GeV resonance as arising from a CP-even Higgs-like particle. 

^It is hard to give a precise unambiguous criterion of electroweak-naturalness, but the general notion is rather clear. 
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show how some of these i?-symmetries are consistent with leptonic or baryonic RPV operators, and 
therefore provide an elegant motivation for the existence of these operators (see 27 for an ahernative 
approach motivating RPV operators). We then focus on the case where the lepton number is tied to 
U{l)ji symmetry, and study it in detail with the goal of laying out the LHC phenomenology. 

Indeed, having an i?-symmetry gives rise to many interesting phenomenological features. For 
example, Majorana gaugino masses, certain scalar trilinear "A-terms", and the "/x-term" are forbidden. 
However, gaugino masses of the Dirac type are allowed]^ This leads to a significant suppression of 
flavor and CP-violating effects relative to the MSSM for 0{1) flavor- violating soft scalar masses and 



phases 29j. Various other aspects of /^-symmetric models have been studied in the literature 30 39 



Although the minimal i?-symmetric spectrum does not give rise to gauge coupling unification, many 



scenarios for adding additional matter have been proposed which could help unify the couplings 28 



40,41 . Particular variants of i?-symmetric models can also give rise to a strong electroweak phase 



transition generating the observed baryon asymmetry, as well as an LSP DM candidate 42 ,43 . Finally, 



a remarkable feature, which we will exploit in this work, is that when the i?-symmetry is identified 
with a lepton or baryon number respectively]^ this allows leptonic or baryonic RPV operators in the 
Lagrangian consistent with these symmetries. Therefore, these are not subject to stringent constraints 
from lepton or baryon number violating observables, such as upper bounds on neutrino masses and 
nucleon-antinucleon oscillations, respectively. 

i?-symmetries are also well-motivated from a more theoretical point of view. In the global su- 
persymmetric limit, it has been known for quite some time that i?-symmetry plays an important 
role in supersymmetry breaking as it is directly related to the existence of supersymmetry breaking 



minima due to the Nelson-Seiberg theorem 44 . Since Majorana gaugino masses necessarily break 
the i?-symmetry, considerable effort has been devoted to generating large-enough Majorana gaugino 
masses while still preserving enough /^-symmetry to keep supersymmetry breaking intact. However, 
an alternative is to consider generic supersymmetry breaking scenarios that give rise to Dirac gaugino 
masses fully consistent with the i?-symmetry. 

The scope of this work is the following. We present a classification of phenomenologically viable 
/^-symmetric models by providing a rather general description of such models in Section |2} This 
is done by making manifest the relevant U{1) symmetries present. Then, we show that preserving 
different combinations of U{1) symmetries gives rise to the different variants of /^-symmetric models 
considered in the literature, including some which have been relatively poorly explored. In particular, 
we will see that in addition to the relatively well-studied i?-symmetric models with the usual "i?-parity 
conserving" operators, there are models which include subsets of the so-called "i?-parity violating 
(RPV)" operators. However, it is important to note that while the standard i?-parity is violated in 
these models, there is still a continuous i?-symmetry and these operators are perfectly compatible 
with it. These kinds of i?-symmetric models arise when one identifies the lepton number (L) or the 
baryon number {B) of the SM fermions with their i?-charges, as will be clear soon. When the lepton 
number behaves as an /^-symmetry, depending upon the region of parameter space, it is possible to 
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Dirac gaugino masses can also be motivated from "supersoft" supersymmetry breaking in which the gauge sector 



has Af — 2 supersymmetry 28 . 

''The lepton and baryon numbers of the SM fermions are standard, but the extension of these to BSM particles is 
"non-standard" . 
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have either the usual down-type Higgs [H^) providing masses to the down-type fermions, or have 
one of the sneutrinos (z/i) providing such masses since it also gets a vacuum-expectation value {vev). 
Furthermore, in the limiting case when Ha is heavy and is not part of the low-energy spectrum, then 
the dominant contribution to down-type masses arises from the sneutrino vev, which is not constrained 
by neutrino masses unlike that in standard RPV models. All these points will be explained in detail 
in the following sections. After describing the general classification of models in Section [2| in the 
rest of the paper we specialize to the case where the lepton number is related to the /^-symmetry, as 
suggested by the title of the paper. We present some characteristic features of the model in Section [3| 
highlighting the differences from standard RPV models. The fermionic electroweak sector of the model 
is studied in Section |4| followed by a discussion of the existing constraints on the model from indirect 
effects, such as electroweak precision observables, flavor physics, etc. in Section |5| The basic aspects 
of the Higgs sector of these models are laid out in Section |6| in particular the region of parameter 
space which gives rise to a ~125 GeV CP-even eigenstate is explained. Finally, we describe the broad 
phenomenological features of the model in Section [7j Since the collider signals of these scenarios 
are quite interesting and novel, in this paper we discuss only the qualitative features which set them 
apart from other models. A more detailed treatment of collider constraints and signals is done in a 
companion paper [l] . Finally, we briefly discuss some aspects of Dark Matter in Section [sj followed by 
conclusions and future directions in Section [9j The appendices deal with some details of i?-symmetry 
breaking operators in Appendix |A| a description of a flavor ansatz for the A and A' couplings (the 
standard notation for the coefficients of the LLE'^ and LQD^ operators, respectively) in Appendix [B| 
and an estimate of the lower bound on A^33 couplings, given the observation of a lepto-quark (LQ) 
signal, in Appendix [C] 



2 Classification of i?- symmetric Models 

We start with the prototypical /^-symmetric Standard Model studied in |29]. In addition to the 
superfield content of the MSSM - Hd, Qi, Uf, D^, Li, Ef {i = 1,2,3), this includes a pair of 
vector-like SU{2)l doublet superfields Ru and Rd (with hypercharge \ and -\ respectively), as well as 
superfields in the adjoint representation of the SM gauge group: a "hypercharge adjoint" or singlet, 
S, an SU{2)l triplet, T (with zero hypercharge), and an SU{'i)c octet, O. 

The relevant global U (1) symmetries of the model in [29] include an i?-symmetry denoted by t/(l)/j(, 
along with the well known lepton number U{1)l and baryon number U{1)b-, as shown below in Tablefl] 
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1/3 


-1/3 


-1/3 






























Table 1: (Super)Field content and U{1) charge assignments. 
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The following superpotential consistent with the above symmetries was considered in 29 : 



It is possible to also write down the following terms with adjoint superfields, consistent with all 
symmetries: 



adj 



S{\HuRd + ^d^uHd) + {X^HuTRd + X^RuTHd) ■ 



(2) 



It is easy to see that with the above i?-charge assignments the usual "RPV" operators, schematically 
denoted by LHu-, LLE'^, LQV^, and U'^D'^D'^, are forbidden. However, /^-symmetries are not incon- 
sistent with subsets of RPV operators in general, since it is possible to construct /^-symmetries Ri 
{i = 1,2,3), which are linear combinations of [/(l)^^, U{1)l and U{1)b, as seen in Table [2] below: 
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f^(l)/?2=-Ro+B 
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U{1)rs=Ro+l 
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1 
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Table 2: (Super)Field content and three different combinations of U{1)r charge assignments. 

We see that depending on the choice of the /^-symmetry, the following RPV operators are allowed 
in the superpotential, in addition to those in ([T]) and ([2]): 



Wi 
W2 



Wo + Wadj + XijkULjEl + X'^j^LiQjDl 
Wo + Wad, + KjiJJtDpi , 



Wo + Wadj+lifH^U . 



{Ri = Rq — L) 
{R2 = Ro + B) 
(i?3 = Ro + L) 



(3) 



In principle, the term ^if^ R^Li is also allowed for R = Ri, but it is possible to do an SU{A) field 
redefinition of the Lj and H^, and define "-ffd" as the field which couples to Ru, leading to 

Thus, this provides a rather general classification of phenomenologically viable i?-symmetric Stan- 
dard Models. The choice R = Ri corresponds to identifying the lepton numbers of the SM fermions 
with (the negative of) their i?-charges, while the choice R = R^ corresponds to identifying them 
with their i?-charges. The choice R = R2 on the other hand identifies the baryon numbers of the SM 



fermions with their i?-charges, and has been considered in 26 . Note that since the i?-symmetries -Ri,2.3 
are identified with lepton or baryon number, they are anomalous and should therefore be thought of as 
accidental low-energy symmetries just like the latter (note that Ro is non-anomalous). The above are 
special limits of a generic U{l)j^ symmetry with R = R + aL + bB for real a and b. However, all other 
cases are severely constrained by proton decay bounds, and we do not consider them any further. 

As we will see later, even though the /^-symmetries Ri,R2 and R3 allow RPV operators, they 
are less constrained than standard RPV models with the same operators. The basic reason for this 
is that in these models these RPV operators are consistent with lepton or baryon numbers (which 



5 



are identified with the i?-symmetries above), hence there are no constraints on the couphngs of these 
operators from processes which violate lepton or baryon number. 



2.1 Super symmetry Breaking 

In order to fully specify the lagrangian, supersymmetry breaking terms must be included. Since we are 
interested in /^-symmetric models, we imagine a situation in which supersymmetry breaking (at least 
in the global limit) is not accompanied by i?-breaking. This can happen if supersymmetry breaking 



is of the D-type, as described in 45 . This includes both D-teim SUSY breaking parametrized by a 



spurion superfield = A'„ + e^D' with R[WJ = 1 {R[X'J = l,R[D'] = 0) and {D') ^ 0, as well as 
F-term supersymmetry breaking parametrized by a spurion superfield X = x + 6'^ Fx with R[X] = 2 
and (x) = 0, (Fx) ^ 0. 

We follow the same procedure as before, by first considering soft terms consistent with the original 
i?-charge assignments (R = Rq) as in Table [T| and then studying additional terms allowed by the other 
choices - {Ri, R2, R3} in Table |2} 

The spurion generates "supersoft" terms via[^ 



d^e — [ciWi^^S + C2W!^^'T + csW^^^'^O"] + h.c. 



(4) 



which contain Dirac gaugino masses iriDi = CiD'/M^. Here denotes the scale of SUSY breaking 
mediation (e.g. the messenger scale in Dirac gauge-mediation scenarios). The above terms preserve a 
U{1)r symmetry under which the and have i?-charge0 1, while R[S] = R[T'] = R[0''] = 0. 

The spurion X generates the following f/(l) ^-preserving renormalizable soft terms: 



- soft 



tsS + ^bsS^ + ^AsS"" + ^brT' + ^boO^ + BfiH^H^ 



+ AsSHuHd + ATHuTHd + A^ST^ + A^SO^ + h.c. 



(5) 



where the sum runs over all the scalars, and we denote the scalar components by the same notation 
used for the superfields. These are generated via the following operators: 



HuHd + e M^S + 52 + + 02 + _L X cubic + h.c 



/ 



X 



d^e j^{ST^ + SO^ + S^) + h.c. 



(6) 
(7) 



We can see that operators quadratic in the visible superfields in the first line in nQv are of order ~^2~- 



^The yvi'^ are the SU{3)c x SU{2)l x U{1)y chiral superfield strengths. 

^The choice i?[yVQ] = RlW'^] = 1 is dictated by the SUSY gauge kinetic terms; in particular this always implies 
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We take Fx ~ D' and 



^ = MsusY ~ 100 GeV - 1 TeV . (8) 

So Dirac gaugino masses from Q, and the non-holomorphic soft mass-squareds, and 5-terms from 
^ are parametrically of the same order (although there may be modest numerical hierarchies). The 
operators which are cubic and higher order in the visible superfields in Eq. ^ will be suppressed by 
powers of Mgusv/^*; and are therefore very suppressed. For the linear term in S in (|5]), dimensional 
analysis generically gives a coefficient ts of order M^M|usy- Phenomenologically however, ts should 
not be larger than Mgygy, since otherwise the scalar singlet tadpole, ts S, will destabilize the hierarchy. 



Ref. 46 has recently argued that this is indeed the case in these scenarios, so that one has e ^ 1. 

The operators in the second line in ^ can give trilinear A-terms involving the adjoint fields of 
order ~ Msusy if allowed. However, these are forbidden if X is not a gauge singlet, implying 
that the scale of these operators can be easily suppressed relative to those in the first line, whose 
scale is naturally set by MgusY- Finally, note that since we imagine that X belongs to a hidden 
sector which has no direct couplings to the observable sector superfields above, there are no terms 
like J X[M^S + Hu Hd + 5*^ + + O^] (due to the non-renormalization theorem, it is technically 
natural to omit these superpotential couplings). 

With other choices for the /2-symmetry - Ri, R2 and -R3, as described in Table [2} one can write 
down additional soft supersymmetry breaking operators as follows: 

^soft _ ^soft ^ BfifH^Li + AfSHuU + A^H^TLi , Ri = Rq - L 

= Lf\ R2 = Ro + B (9) 

T soft r soft p D I r 

L,^ — Lq . Ks — Kq + L 

Only the case Ri = Rq — L allows additional gauge-invariant operators consistent with the i?-symmetry 
(since L has zero i?-charge). These additional soft terms for the case i? = i?i in (g give rise to a very 
interesting possibility for the vacuum structure of the theory. 

The presence of the -B/x^*'' term in (9) for R = Ri implies that one of the left-handed sneutrinos gets 
a fcf [^due to a tadpole for the sneutrino when Hu gets a vev. It then becomes possible to distinguish 
two extreme cases: i) (z>i) ^ (H^), and ii) (z>i) 3> (H^)- In fact the size of the fidRuHd term in 
superpotential Wq controls which one is relevant. This is because, schematically. 



Here mi is the soft mass of the left-handed sleptons. Thus, if /i^ ^ m|, then (z/i) ^ {H^)- 

^Recall that we have defined as the hnear combination of SU{2)l doublets with Y = —1/2 and R = 0, that 

couples to Ru in the superpotential. Within this class of bases, it is further possible, by SU{3) rotations, to go to the 
"single- vev-basis" where only one of the three Li acquires a vev (see Section [s]) . We will see later that the most natural 
choice is to identify the direction of this vev with the "electron" direction, i.e. i = 1. The Hd vev in this basis may be 
non- vanishing, but we are interested in a region of parameter space where it is small compared to the EW scale. 
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In the remainder of this paper, we will study the choice R = Ri and the case (z/i) ^ (H^) in detail 
as it gives rise to rather novel and interesting phenomenology. We will discuss the phenomenology of 
the R = R2 case in Section 7.6 very briefly, since that case has already been considered in ^6]. For the 
case we are interested in, the superpotential in ([T]), (|2]) and the soft terms in (|5]), ([9]) can be simplified 
since the large fid term allows us to integrate out the fields Ru and H^, leading to the following: 



W 



'soft 



l^uHuRd + ^u^HuRd + ^u^uTRd + y'j'jHuQiUj + XijkLiLjE'^ + X[ji.LiQjDl 



tsS + ^bsS'' + hrT^ + hoO' + Bfif H^U 

+ \AsS^ + A^tST"" + A^oSO^ + AfSHuLi + A^HuTU + h.c. 



:iii 



:i2i 



where the terms in the second line in Lgoft are assumed to be suppressed relative to MgusY for simplicity, 
from the arguments below ([T]). Note that in the R = Ri scenario the down-type masses arise from the 
LLE'^ and LQD'^ operators when the left-handed sneutrino gets a vev (assuming (H^) ^ (^'i))- 



2.2 i?-breaking 

It is well-known that the vanishingly small value of the cosmological constant breaks i?-symmetry since 
it requires a non-zero value of the superpotential in the vacuum, and the superpotential has non-zero 
i?-charge. Since the gravitino mass ~ (W) (in Planck units), this implies that m3/2 is the order 
parameter of i?-breaking. As mentioned in the Introduction, in this work we imagine a setup in which 
m3/2 is much smaller than the TeV scale. Hence, the effects of i?-breaking are also small. 

The breaking of i?-symmetry will eventually be transmitted to the visible sector. This can essen- 
tially happen in two ways. A simple possibility is that i?- breaking is mediated to the visible sector by 
generic Planck suppressed operators, which we denote as "generic gravity mediation" . This is a natu- 
ral possibility since gravity is expected to violate all global symmetries in general. However, another 
possibility is that these generic Planck suppressed operators respect the /^-symmetry (at least to a 
very good approximation) due to it being an accidental symmetry of the visible sector, see [47j for an 
example. In this case R breaking will generically be communicated to the visible sector via anomaly 
mediation. 

In fact, the source of i?-breaking can be connected to observable physics in an interesting way. For 



example, for R = Ri the breaking of the i?-symmetry will give rise to neutrino masses 24 ,25 , while 



for R = R2, i?-breaking will give rise to nucleon-antinucleon oscillations and may also lead to proton 



decay in certain cases 26 . Existing constraints from these observables then put an upper bound on 



7713/2 and therefore on the messenger scale, M^:, as well 24 . In Appendix 2.2, we describe some details 
of the sizes of i?-breaking operators. However, the collider phenomenology is largely determined by 
the approximate i?-symmetry, and therefore we will often focus on the /^-symmetric limit. A more 
thorough analysis of the full effects of /^-breaking is left for the future. 
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3 Characteristic Features of the Model with R = Ri 



In the i?-symmetric limit, the superpotential and soft terms are given by Eqs. (11) and (12). The 



terms including the lepton (super) fields above are written in a general basis. In a general basis, all 
the sneutrino fields can develop non- vanishing vacuum expectation values (vev^s) since there are Bfi^ 
terms for all of them. However, it is always possible to choose a basis in which only one of the sneutrino 
fields gets a vev 25 ,^ In this sense, there is a similarity with models of i?-parity violation (RPV) where 
the Hd and Li fields mix, in general, with each other, giving rise to mass eigenstates La, a = 1, ...4. 
There is an important difference, however, since has been integrated out, and the "light down-type 
doublet", Rd, has a different i?-charge from the Lj. Thus, there is only a three-dimensional space in 



which the (sneutrino) fields can mix. We will follow the analysis in 21 keeping this difference in mind. 
The general basis can be related to the "single vev basis" as: 



Via) 



L(a) + ^ eibLh . 



(13) 



Here the lepton of fiavor (a) is assumed to get a vev, while k run over all three generations and b 
(and later c) runs over only two generations (those which do not get a vev). The ctb are the matrix 
elements that relate the fields in the two bases. There is still the freedom to rotate Lb, and by choosing 
an appropriate Cib one can go to a basis in which the charged Yukawa couplings are diagonal. 

Since the lepton Yukawa coupling is provided by the Xij^LiLjE^ operator, and Xijk is antisymmetric 
in the first two indices, gauge invariance prevents the lepton of fiavor (a) from getting a mass from such 
operators. Its mass can, nevertheless, be generated from supersymmetry breaking (but /^-preserving) 
operators. For example, it could come from the following operators: 



^(1) jd'0^,HiL,^,El,, 



II) 



y{2) 



c 

(a) 



M2 



Am 



(1) 



Vii) 



x\ 



SUSY 



Am 



(2) 



, \Fx? Vu 
^(^) Ml - V 4/.2 



(14) 



We see that the first operator can provide viable lepton masses only if the messenger scale is 



low 124 , 48 , but the second operator is generated after integrating out the fields Ru and Hd with a 



large supersymmetric mass term HdRuHd (see discussion around (10), and appendix B in ^2]), and can 



give a viable contribution even if the messenger scale is high. Note that both these operators are 
present for all lepton fiavors in general, so these will provide contributions to lepton masses in addition 
to those from the superpotential in ( 11 )|^ Therefore, the smallness of the electron mass makes it natural 
to take (a) = 1(e), and 6, c = 2 (/i), 3 (r), so that the electron gets its mass solely from supersymmetry 



breaking (but i?-preserving) operators in (14) 25 . We will assume this henceforth. 



^The physics is of course basis independent. 

^They will also contain flavor ofi^-diagonal entries, presumably of order me 



In the lepton mass eigenbasis, these 
will induce off-diagonal slepton masses, even if these are diagonal (but not degenerate) in the gauge eigenbasis. In this 
i?-symmetric framework, we expect the constraints from /i ^ 57 and fi — e conversion in nuclei to be satisfied since the 
mixing angles are of order me/m^ and because we take the Dirac bino mass Aff around 1 TeV (see section [?]) 49 . 
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In the "single-vev" and "mass-eigenstate" basis, the superpotential is given by: 



W — HuHuRd + ^u^HuRd + ^u^uTRd + WVukawa + Wrrilinear , 

M^Yukawa = 2/i^^L(,)L,E,^ + l/f^(a)4A^ (15) 

b=2,3 i=l,2,3 

W^Trilinear = A23iI/2-f'3-E'i^ + Kij^bQiDj . 

j=l,2,3 ij=l,2,3;fe=2,3 

where the hat denotes that the lepton and quark fields are in the "mass-eigenstate" basis. Note that 



the first two indices in the trihnear term LLE'^ in PVrriiinear in Eq- (15) are fixed to be (2) and (3) 
since (a) = (1), and since the couphng is antisymmetric in the first two indices. There is, however, no 
such antisymmetry for the first two indices in the LQD'^ term in M/rriUnear- Overall, these terms have 



a rather different fiavor structure compared to analogous trilinear RPV couplings in RPV models 21 

To contrast some other important features of our model against RPV models considered in the 
literature, it is instructive to look at properties of standard RPV models with both bilinear and trilinear 



RPV operators, where we use the established results summarized in 21 . For example, the presence 
of the fermionic and scalar bilinear RPV operators in such models is associated with mixing between 
neutrinos and neutralinos (or charged leptons and charginos). In particular, the superpotential bilinear 
/i^ HuLi generates neutrino masses at tree-level proportional to tan^ ^ through such mixings, where the 
angle ^ parametrizes the physical higgsino-lepton mixing in the fermion sector (which cannot be rotated 
away) in a basis-independent way. This leads to a very stringent bound, sin ^ < 3 X 10- Vl + tan2 /3. 
Hence, in the basis with a single /i-term (i.e. fiHuHd, but no ^'^HuLi), the sneutrino vevs are forced 
to be extremely small due to the upper bound on neutrino masses. Furthermore, the presence of the 
scalar bihnear RPV operator Bix^^HuLi in Csoft and the trilinear RPV operators XLLE'^, X'LQD'^ in 
the superpotential, give rise to one-loop contributions to neutrino masses proportional to i) g'^ {Bfi^l'y, 
ii) gX' (or gX), iii) A'^ or (A^), which put stringent bounds on many A, A' couplings (such as Ajsa 
and A^33, i = 1,2,3) as well as the size of B^l (this can also be interpreted as putting a stringent 
bound on the sneutrino vev^s through the basis-independent Higgs-slepton mixing angle sin^, even in 
the absence of the superpotential bilinear RPV operator). 

In the model under consideration, the operators XLLE'^ and X'LQD'^ in the superpotential, and 
the operator BjjifjHuLi in Csoft, preserve a lepton number (which is identified with the i?-symmetry 
R = Ri), unlike that in standard RPV models above. Hence these terms cannot generate Majorana 
neutrino masses which violate the lepton number {R = Ri) either at tree-level or loop-level, as long as 
the i?-symmetry is conserved. This further implies that the sneutrino vev (induced by the -B/i^'' term) 
can be significant and can play the role of a Higgs field. Also note that there is no fi'^H^Li term in 
the Lagrangian, implying that the basis in which the Yukawa couplings in VTvukawa are diagonal is the 



same as the mass-eigenstate basis of the charged leptons, unlike in RPV models 21 



Finally, from above we see that the bounds on neutrino masses are only relevant when /^-breaking 



effects are taken into account and are, hence, proportional to 171^/2 (see Section 2.2). Thus, for given 



values of the sneutrino vev and the A, A' couplings (which are consistent with other constraints, see 
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Here g schematically denotes any of the two electroweak gauge couplings. 
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Section [5|, the bounds from neutrino masses only provide a bound on 7713/2- 

4 The Fermionic Electroweak Sector 

Since in our framework lepton number is identified with an /^-symmetry {R = Ri), the neutrahnos and 
neutrinos on the one hand, and charginos and charged leptons on the other, share the same quantum 
numbers (in particular, their i?-charge). They can, therefore, mix after electroweak symmetry breaking. 



4.1 Charginos & Charged Leptons 

The charginos and charged leptons will mix after electroweak symmetry breaking in general. However, 



from (15), it is clear that only the charged lepton of flavor (a) (the electron since (a) = 1) will mix 
with the charginos since only that flavor gets a vev. One has four Dirac charginos (one of which is the 
electron), which can be further split according to their electric and i?-charges. Then, one can form 
two groups of 2-component fields - one with R = +Q, and one with R = —Q. This implies that the 
chargino mass matrix can be written as: 

= {w^,e%), r- = (T,-,eZ), (fori? = + Q) 

r+ = (w-R,), = (T+,iJ+), (fori? =-Q) 

where 

M?' = . M^,:' = ^l'-"' ) . (17) 

^ \ me J ^ \gvu -flu- X^Vs + XIvtJ 

The notation x~^~ , for instance, implies that the field x has electric charge +1 and i?-charge -1. Here 
stands for the Dirac wino mass, g is the SU{2) gauge coupling, me is the electron mass, and Vu, 
Vs & vt are the vev's of if°, z>(a), 5" and T° respectively. The above matrices can be diagonalized by 
two pairs of 2 x 2 matrices - {V"*", U"*"} for (x"*"*" X )> ^'^^ {y~ ) U~} for {x^~ x~^)i respectively, such 
that (V''")"I'M[^^U"'" and (V^)"I'M|-7''U^ are diagonal (and with positive eigenvalues). The above states 
are naturally arranged into four 4-component Dirac fields X^^ = {xf^Xi ) and Xl^" = {xt^xi^), 
with i = 1,2, whose charge conjugates are denoted by X^ and X^^^ respectively. In this notation, 
e = X'^ corresponds to the physical (Dirac) electron field. 

4.2 Neutrahnos & Neutrinos 

Similar to the charged fermion sector, the neutrahnos and neutrinos in the neutral fermion sector will 
mix after electroweak symmetry breaking. For simplicity, we consider only one neutrino generation 



since adding the other neutrinos does not change the qualitative collider picture (see 25 for a thorough 
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discussion of neutrino masses and mixings). Also, since Majorana neutrino masses violate lepton 
number (hence R symmetry in our case), they are massless in the /^-symmetric limit It is convenient 
to write the matrix in the "Dirac" basis by grouping the fields with i?-charges +1 and -1 separately. 

Then, similar to the charginos, the mass terms for the neutralinos can be written a.s Cn = 
(X°+)^ where = (6°, R^^) and x°" = {S, T°, H^, z/e). The notation implies that the 

field X has vanishing electric charge and -R-charge +1. The mass matrix M^r is given by: 



M 



N 



D 



I Mi 




a/2 
_gvu 



V2 

gV(a) 





;i8) 



After diagonalizing the above mass matrix by unitary transformations and U^, as in the chargino 
case above, one obtains three Dirac mass eigenstates = {x^i^ X^i')^ with i = 1,2,3, and one 
massless Weyl neutralino x^~ that necessarily remains massless, which is identified with the massless 
neutrino eigenstate, and is given in general by: 



ug~s + u!:if' + uzm + u^ 



N . 



it 



(19) 



By a slight abuse of notation, sometimes we will also refer to Xa "^'e", where it will always refer to 
the mass eigenstate, and should not be confused with the original gauge eigenstate. 



5 Indirect Constraints 



The identification of the lepton numbers of SM fermions with (the negative of) their i?-charges is 
subject to various "indirect" constraints. As explained in Section [3} the sneutrino vev can be significant 
in these models since there are no bounds from neutrino masses. Thus, the most stringent constraints in 



these models arise from electroweak precision measurements and from fiavor violating processes 24,25 



We will see that these can be satisfied in a large range of parameter space and for reasonable flavor-off 



diagonal couplings. Many of the constraints are similar to those studied in 24 . It is convenient to 
divide the constraints into two categories: 

• Constraints on the sneutrino vev or equivalently tan/3 = 

• Constraints on the A and A' couplings. 



5.1 Constraints on the sneutrino vev V(^a) (tan^) 

First, the mixing between the charged leptons and charginos gives rise to a deviation in the couplings 
of the Z to charged leptons in general, which is constrained by electroweak precision measurements. 
As shown in the previous section, since only one lepton (of flavor (a)=l) mixes with the charginos, the 

"'^^We are assuming that there are no right-handed neutrinos with i?-charges such as to aUow writing down Dirac 
neutrino mass terms consistent with the i?-symmetry. 
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mixing is identical to that in 24 , giving rise to the following deviation in the vector and axial- vector 
coupling to the Z from those in the SM: 



sm 



(20) 



sm I 



Aml{M^Y 



2gv^a)M, 



D 



From the measured value of the coupling g\ = —0.50111 ± 0.00035 [50j, one gets an upper bound 



on I TTu )• For example, for M^= 1500 GeV (300 GeV), the current data puts an upper bound 
"^(a) ^ 61 GeV (12 GeV) at la. In this work, we assume that the (Dirac) gauginos are heavier than 
the scalars. As a benchmark, we take ^1.5 TeV henceforth, implying that v^^^ ~ 60 GeV. The 
fact that only the charged lepton of flavor (a) (e in our case) mixes with the charginos also gives rise 
to constraints from charged current universality. However, the bounds from these are not as strong as 
those derived from the Z-coupling above (see [24]). A lower bound on f however, arises from leptonic 
Yukawa couplings, t/r in particular. This is because the leptonic Yukawa couplings arise from the LLE'^ 



operator, and hence are part of the A couplings [see Eq. (11)]. Therefore, these give rise to extra tree- 
level contributions to electroweak observables similar to those in traditional RPV models (for a review 
of constraints in RPV models, see (21]). One finds that the dominant constraint arises from the r 
Yukawa coupling [y^. = y^^'' = A133) contributing to the ratio = T{t — )■ e z/g /T{t ^ ji 
as shown in Fig. [l] (i?^ is normalized to the dominant r decay due to W exchange). This gives: 



24 



yr < 0.07 

which puts a lower bound V(a) ^ 8 GeV (2.5 GeV) 
for mfj^ = 300 GeV (1 TeV). Since in our frame- 
work, one of the sneutrinos behaves as a Higgs field 
which is expected to be around the electroweak 
scale, it is natural to expect that the masses of all 
sleptons are around the electroweak scale, i.e. few 
hundred GeV (see more discussion on this in Sec- 
tion [7j). We, therefore, take v^^"- ~ 10 GeV. Com- 
bining with the upper bound on V(a) from above, we 
obtain the range: 



100 GeV/ 



(21) 



T 



L 




y. 



T 



Figure 1: Contribution to the decay width ^■ 
eJV^VT, where both the interaction vertices correspond 
to the r Yukawa couphng Ht- = A133. The arrows indicate 
i?-nuinber flow. 



10 GeV < v^a) < 60 GeV 
or 17.4 > tan/3 > 2.7, 



(22) 



where tan/3 = fu/f(a), which we use in our subsequent analysis. Thus, we see that the sneutrino vev 
can be much larger than in standard (bilinear) RPV models. This is one of the most distinctive features 
of the model, and it also plays a crucial role in LHC phenomenology as we will see. In Section [7| we 
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will discuss how the existence of such a large sneutrino vev could be inferred at the LHC 



5.2 Constraints on A and A' couplings 
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As explained at the end of Section |3} in our model the bounds from neutrino masses can be interpreted 
as a bound on the gravitino mass (since it is the order parameter of -R-breaking) , implying that the 
bounds on A, A' couplings only arise from other observables, such as flavor- violating processes. We 
will see that this has an important effect on the bounds on the A^ss and A'j33 {i = 1, 2, 3) couplings in 

^33 < 10-3) on these 

On the other hand, the upper bounds on Ai33 and A^33 in our model arise 
from flavor-violating processes and can be rather mild as we will see. 

Constraints from flavor-violating processes in the lepton and hadron sector provide bounds on the 
A and A' couplings. Although there do exist bounds on single couplings, they are typically weak, and 
most of the stringent bounds arise from the products of two couplings. These can be classified into 
three categories, those constraining i) A A couplings, ii) A A' couplings, and iii) A' A' couplings. We 



particular, because in standard RPV models the most stringent bounds (Aj33,A^ 
arise from neutrino masses 



take our results from the general analysis of constraints arising from RPV operators in 21,51 53 
However, as mentioned above, in our framework the lepton and down-type Yukawa couplings are part 
of the A and A' couplings respectively, and the results of 21,51 53 must be interpreted accordingly. 
In particular, using the notation in (15), A122 = 1/2'^ and A133 = |/^% while A'^^^^ = y['^\ X'122 = y2^\ and 
A']^33 = y^/K The remaining A, A' couplings are included in PVrriiinear in (15). Among these non- Yukawa 

= 1, 2, 3), and 
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- A 



23i 



trilinear couplings, there are only three independent couplings of the A-type 
two 3x3 matrices of the A'-type - A2j^ and X'^j/^ with independent entries. Thus, not only the number 
of independent couplings is greatly reduced, but also the implications of these bounds on the parameter 
space of the model are very different compared to those within a standard RPV scenario. 



Coupling(s) 


Upper Bound 


Process 




23fc 


;k=l,2,3 


0.07 CfcR 


[Rr, Rtii] 21 


51 






A231 Z/a 
or, IA231 




5.2 X 10-4 [i>L,f 
0.052 cos 13 [i)L,f 


[r — 7- e e /i] 


52 






\ (e) 
^232 Us 

or, A232 




7.0 X 10-4 [i)L,]' 
0.070 cos (3 [i>L2? 


[r — )• /i e /i] 


52 






X (e) 
A233 ^3 

or, A233 




2.2 X 10-4 
0.022 cos 13 


[t — eP°/ /X — e in nuclei] 53 



Table 3: Some upper bounds on A couplings from various flavor-violating processes. Here [i^La]^ = ( looGcv )^- "^^^ 
notation for the processes is the same as in the corresponding references. 



^^In standard RPV models, the contributions from A, A' couplings are proportional to the quark and lepton masses, 
hence the bounds are rather tight for the third generation (s)quarks: A^sa and A-gg, for i = 1,2, 3. 
""^^Note that the first two indices in the A couplings are anti-symmetric. 
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Coupling 




Upper Bound(s) 


Process 




A211 V2 ^ 
or, IA211 




2.1 X 10-s [z/L,]2 
3.5 X 10-5 cos/3[z>L2]^ 


[/i — e in nuclei] 


21 


,51 






A21 
or. 


22/2 
\/ 

^^212 




6.0 X 10-9 [i>L2]2 
1.01 X 10-5 cos/3[z>L2]2 


[f^^ /ie/e/i] £ 


2 




A213 V2 '* 
or, IA'213 




1.3 X 10-5 \vi^^^ 
0.022 cos/3[z>L2]^ 




52 






A22 
or. 


(e) 
l2/2 

\/ 

^^221 




6. X 10-9 [z>L,]2 
1.01 X 10-5 cos/3[z>i2]^ 


[i^o ^ /ie/e/i] |52] 




A' v^^^ 

^222 2/2 
or, A222 




1.0 X 10-5 
0.032 cos/3 


[r — i- eP^/ /i — e in nuclei] 


53 






^^223 2/2 

or, A223 




7.6 X 10-5 [i>^^]2 
0.128 cos/3[z>i2]^ 


[5° ^ e/i] 


52 




-^231 2/2 ^ 

or, 


! A231 ?/3 ^ 
A23I 


1.3 X 10-5 \vl-^^] 1.6 X 10-3 [ml3]2 
0.022 cos/3[z>L,]2; 0.099 cos/3[ml3]2 




52 




^232 2/2 

or. 


1 1 ^^232 2/3 1 

X' 1 

^232 1 


7.6 X 10-5 [i/La] , 2.7 X 10-^ [ul.^]^ 
0.128 cos/3[z>L,]2; 0.016 cos /3[uls? 


[5° ^ /ie] 52 ; [fe s/ie] 


51 






\/ 

"^23 

or. 


(d) 
32/3 

\/ 

^233 




1.1 X 10-5 
6.8 X 10-3 cos/3 


[r — i- eP'^/ /i — e in nuclei] 


53 





Table 4: Some upper bounds on ALj, (j, k — 1,2, 3) couplings from various flavor-violating processes. 



Starting with the A couplings, we find the bounds in Table [3j We have only listed single coupling 
bounds, and those product (A A) bounds in which one of the couplings is a Yukawa coupling. This 
is because the Yukawa couplings are known up to tan /3, hence they provide the most robust bounds 
(however, we do consider the full set of constraints in our analysis). The A' couplings are more 
numerous. Again, upper bounds exist on various products of the form A A' and A' A', and are listed 



in |21 ,51 53 . As before, we describe the bounds on those products, one of which is a Yukawa coupling 
(either of the leptonic type or the down- type quark type). For some couplings, bounds exist from more 
than one experiment. We list the dominant bound in such cases, unless the bounds are comparable in 
which case we list all of them. The bounds for Ag^^ and Ag^j^. are listed in Tables |4] and [S] respectively. 
As explained earher, the bounds on Ajss and A^33 {i = 1,2,3) are relaxed relative to those in standard 
RPV models. In particular, the bound on A333 is very mild from Table it can be comparable to the 
electroweak gauge couplings. This can have an important effect on collider phenomenology related to 
the third generation (see Section [T]). 

In order to get a better idea about the constraints on these couplings, it is useful to understand what 
generic expectations we have for the spectrum of the model. This will be discussed more in Section [7j 
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Coupling 




Upper Bound(s) 


Process 




or, IA'311 




8.5 X 10-^ [Dl,]'^ 
8.47 X 10-3 cos/3[z>L3]2 


[t — e?]] 


52 






A31 
or. 


22/3 
\/ 

^^312 




9.7 X 10-4 [Hf 
0.097 cos/3[z>L3]2 


L ''J 


52 






A313 yi ^ 

or, IA'313 




3.7 X 10-4 [z>L3]2 
0.037 cos/3[z>L3]2 


[B', ^ ef ] 


52 






A321 Vs ^ 
or, A321 




9.7 X 10-4 [uLf 
0.097 cos 13 [ul;? 




52 






A322 yi ^ 

or, A322 




4.6 X 10-4 [i>^j2 
0.046 cosP[uls? 


[r ^ er^] [52j 


or, 


'1 A331 1/3 
A33I 


3.7 X 10-4 [z/La]^; 2.7 x 10-^ [MLg]^ 
0.037 cos/3 [i>L3]2; 0.168 cos/3 [ml3]2 


[B', ^ re] 


52 






^333 i/3 
or, A333 




2.1 X 10-2 
1.305 cos/3 




53 





Table 5: Some upper bounds on Ag^-j, (j, k = 1,2, 3) couplings from various flavor-violating processes. We do not show 
bounds on A323 and A332 as they do not appear in product bounds where the other coupling is a Yukawa coupling. 



here we just make some brief remarks. We imagine a situation in which the (Dirac) gauginos are 
heavy and the scalars and Higgsinos are relatively light. For concreteness, we take /i = 200 GeV, 
m| ~ m| ~ (200 - 300GeV)2,Mf = 1000 GeV, ~ 1500 GeV. The masses of the squarks are 
determined from current LHC constraints, which are studied in detail in the companion paper [T]. It 
turns out that the bounds on masses of the first two generation squarks are in the 600-700 GeV range, 
while the bounds on the third generation squark masses are lower, around 400 GeV. 

With the above spectrum in mind, it is straightforward (although tedious) to check that (almost) 



all of the remaining bounds on the products of A, A' couplings in |2T|51 53] can be satisfied if the 
values of the A, A' couplings are assumed to saturate the bounds in Tables |3| |4] and |5| The exception 
is the following product bound: 



Coupling(s) 


Upper Bound 


Process 


IA231 A31J 


2.1 X 10-8 [i>L,? 


[r — > e P^/ — e in nuclei] [ss] 



A simple choice, therefore, is to assume that the coupling X'^^^ is negligible (^ 0), while all the other 
couplings still saturate the bounds in Tables |3||4] and [5j We will see in Section[7]that this leads to rather 
interesting prospects for various signals such as lepto-quark-like signals, and single slepton/sneutrino 
production. Finally, it is worth mentioning that it is possible to make other (reasonable) ansatze about 
the flavor dependence of these couplings. In Appendix [B| we discuss another simple ansatz about the 
flavor dependence of these couplings which allows a nice understanding of the relative magnitudes 
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of the various couplings and also satisfies existing constraints. It also has the advantage of further 
reducing the number of independent A and A' couplings. 



6 The Scalar Electroweak Sector 



In this section, we discuss some important features of the electroweak scalar sector of the model, and 
compare and contrast it with the well known case of the MSSM. This is of considerable importance, 
especially after the recent discovery of a Higgs-like particle near 125 GeV (2||3]. A more thorough 
treatment of these issues, including the couplings of the Higgs-like particles in the model, will be the 



subject of another work 54 



From Eqs. (11) and (12), the electroweak scalar potential takes the form: 



EW 



irEW 



V 



EW 



m 
1 



Vp V £) -t- Vsnft "T V'l, 



dW 



soft ^ loop ' 
3 



2 

Y 5 



(23) 
(24) 



a=l 



^ 1 1 



while ^oop refers to the one-loop contribution to the effective potential, which will be specified below. 
Here i runs over all the electroweak scalar fields which could receive vev^s, and Dy and are the 
hypercharge and SU{2)i D-terms, respectively. Compared to the MSSM case, the D-terms contain 
additional pieces associated with the SU{2)l and U{1)y adjoint fields: 



7-ja 



D 



Y 



g{HyHu + R^Rd + lyii + TU'^T) + (Mf T"" + h.c.) 
^-^{HlH^ - R^R, - L\U) + V2 {M^S + h.c.) , 



(25) 



where r'* and A'' are the two and three-dimensional SU (2) generators respectively. Note that the D 
terms above give rise to new trilinear couplings in the scalar potential. Also, the masses of the real 



and imaginary parts oiS = SR + iSi and T = Tr + iTj are split in Eq. (23). For instance, if Mf, bs 



and are real, then = + bs + 4(Mf')^ and m| 



Si 



and rri^^ 



— bs while rriip^ 



ml + bT + 4(Mf )2 



— bx- For simplicity, we assume that there are no CP- violating phases in the potential. 



In order to minimize the above potential, we point out some important simplifications. First, EW 



precision constraints on the p-parameter require the triplet Higgs vev, (T^) 



to be small (< 3 



GeV [55]), which is naturally achieved if the triplet soft breaking mass ^ TeV. Therefore, the 
effect of the triplet on the minimization of the potential must be small, and vt can be set to zero in 
the first approximation. Second, since the i?-symmetry forbids the term B^HuRd, it is easy to see 
that (Rd) = if m^^ > 0, i.e. there is no spontaneous breaking of the U{1)r symmetry. Also, because 
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Rd has a different i?-charge (= 2) than the rest of the electroweak fields (= 0), the degrees of freedom 
in Rd do not mix with those in the other fields and decouple from the rest. 

It is important to understand the similarities and differences in the structure of the scalar potential 



relative to that in "supersoft" SUSY breaking studied in 28 . Since the gauge sector of the model is 



similar to that in 28 , the model shares the good feature that unlike the MSSM, the usual logarithmic 
divergence from the stop contributions to the Higgs mass-squared parameter m|^^ are cutoff by the 
Dirac gluino mass, leading to only a finite contribution. Thus, in contrast to the MSSM, a Dirac gluino 
mass in the multi-TeV range is consistent with electroweak-naturalness. 

On the other hand, the matter sector of the model is rather different from that in [28 1. Indeed, in 
the latter case one gets a vanishing D-term contribution to the Higgs quartic coupling at tree level, 
which is obviously not a good starting point to obtain a Higgs mass near 125 GeV. In our model, 
however, the above conclusion is circumvented by the presence of soft (but not supersoft) operators 
arising from F-terms of the X spurion in ^ which yield the soft parameters {771^,771^,65,65"}, etc. 



Also, the presence of the superpotential couplings in (11) proportional to and A„ can give rise to 
new F-term contributions to the Higgs quartic coupling at tree level if Rd gets a vev, as in 56 , 57 . 
However, since in our model (Rd) = 0, this tree-level F-term contribution is not present. Nevertheless, 
the Af and couplings do provide important contributions to the Higgs quartic coupling at loop 
level. We will see in the next subsection that this is very important in obtaining a CP-even mass 
eigenstate with mass ~ 125 GeV. 

With the above simplifications, it suffices to minimize the scalar potential with respect to the 
neutral fields - t'(a), Sr} to study electroweak symmetry breaking (EWSB)|^ 

^ /x^ + - l^ ^'^^y^'^ ^ v\,p + (gMAa) ^2^2 ^ ^g'^^M^ + Ag7;.(2/i + Aft;,) + r^'B^^f, 



m 



^ + (Af )^ vH}] V, - ^^Mfv'c2, + {ts + Af /i v'sD . 

Here, Sjb stands for sin/3 and so on, and {6Xu, 6X^,6X3} denote the dominant radiative corrections to 
the quartic terms: ^6 Xu {\H'^\'^y , |5Aj, (|i>(a)p)^ and ^(5A3 |i/°p|z/(a)p respectively. In the limit where 



A^ is negligible, the leading-logarithm contributions to these radiative corrections are given by 56 



sx^ ^ M log ( "^h \ ^ ^JgT f ^ 

167r2 \ ml J 167?^ y v'^ 
^ M log to) + ^ log (4 I . (27) 



where the renormalization scale is taken to be close to the electroweak vev. The trilinear soft terms 
A\ have been neglected here since they can be suppressed for reasons mentioned below (pi). In the 



*For numerical results we do a full analysis, including all uew's, and based on the full Coleman- Weinberg potential. 
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above approximation, the CP-even and the CP-odd neutral Higgs fields are linear combinations of the 
real and imaginary parts of {H^, i'l, S} respectively. The charged Higgs H^, on the other hand, is a 
combination of {if^,e|^}. 



6.1 The ~125 GeV Eigenstate 



It is important to understand what region of parameter space of the model gives rise to an eigenstate 
with mass near 125 GeV, given the recent discovery of a Higgs-like particle with that mass. We will 
only make some general and somewhat qualitative comments, leaving a detailed study of these issues 



for future work 54 



To start, let us write down the tree level {6Xu 

Sji) basis: 



5Ay = 6X3 = 0) mass matrix for the CP-even 



neutral states in the 



H 



9' 



,W1 



\vsp[V2g'Mr + 2Af (Aft;, + 



2^Bf,[^^] 

D 
1 



vsp[V2g'Mf + 2Af(/i, + X^Vs)]\ 
-V2g'vc^Mf 



2vs 



where we have used the minimization conditions in (27) to get rid of the non-holomorphic soft mass- 
squareds for U(^a) and S. 

By inspection, one can see that 
the mixing angle between Sr and 
{if°,z/(a)} is essentially controlled 
by the ratio Vs/v. Hence, a larger Vs 
will make this mixing angle larger, 
pushing down the lightest eigen- 
value due to "eigenvalue-repulsion" . 
Thus, Vs should be small in order 
to maximize the lightest eigenvalue 
(we do not consider the possibility of 
a very light singlet scalar). In this 
limit where the off-diagonal entries 
are relatively small, it is then not 
hard to see that the largest eigen- 
value is predominantly Sr, while the 

— z>(a) block gives rise to a tree- 
level smallest eigenvalue approach- 
ing that in the MSSM. 

Thus, in order to obtain the lightest CP-even Higgs mass around 125 GeV, a reasonably large 
radiative contribution to the Higgs quartics (primarily 5A„) is required]^ However, unlike the MSSM, 



Benchmark I 


Benchmark II 


tan /3 = 3 


tan/3 = 17 


Af = 0.1 

>Z = 1-0 
/x„ = 200 GeV 
= 200 GeV 
Mf = 1000 GeV 
BfiL ^ -(174GeV)2 
ts ^ (174GeV)3 
~ (1115GeV)2 
~ (1450 GeV)2 
m|=m| = (500 GeV) 2 


Af = 0.1 
A^ = 0.9 

= 200 GeV 
Mf = 200 GeV 
Mf = 1000 GeV 
BfiL - -(123GeV)2 
ts - (138GeV)3 
~ (880 GeV) 2 
ml ~ (1390 GeV)2 
ml = m| = (500 GeV)^ 


rrih ~ 125 GeV 


rrih ~ 125 GeV 



Table 6: Two benchmarks giving rise to a lightest CP-even Higgs 
mass close to 125 GeV. We take = —5 GeV in both cases. 



""^^As explained earlier, in this model there are no additional tree level contributions proportional to (Af )^ or (Aj)^, 
unlike that in 56 57 . 
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where the dominant contribution to SXu is provided by the stop squarks, here the adjoints S and T 
can also provide a significant contribution through terms proportional to and in (11 ). In fact, it 



is possible that the bulk of the radiative contribution is provided by the triplets T in the loop, with Af 
small and A^ close to unity [see the approximate expressions for {SXu, SX^, 6X3} in (27)]. We give two 



benchmark examples in Tabled with values of the important parameters and for two choices of tan /3, 



close to the minimum and maximum values in (22). In these benchmark examples, the stop squarks 



are taken to be around 500 GeV, so the dominant radiative correction is provided by the triplet T. 
It turns out that even though the singlet and triplet scalar masses are > TeV, the sensitivity of the 
Higgs potential on them is less than that on heavy stops in the MSSM which could generate a ~125 
GeV Higgs mass, making this model significantly less fine-tuned than the MSSM. Part of the reason 



is to be found in the factor of 5 versus 3 displayed in (27), which makes the contribution to the Higgs 
quartic more effective for the triplets than for the stops. For instance, for yt = A^, a 500 GeV (1 TeV) 
triplet gives the same contribution to SXu as stops with m^^ = m^^ = 1 TeV (m^^ = m^^ = 3.2 TeV). 
In addition, one finds that the radiative correction of the triplet to m'jj^ corresponds to that of a single 
stop with m? = and yt = A?^. Hence, for fixed SXu and equal stop masses (no stop LR mixing), 
one can estimate an overall suppression of the triplet "quadratic divergence" compared to that of the 
stops by a factor of about 



^^Jtripiet 1 (AD' e^'^"'* 



Till of n-r\o '-' 



{2i 



I stops ^* g 3b| 



For A^ ~ 1 and tjit ~ 1 TeV (corresponding to SXu ~ 0.11), the above represents a suppression by a 
factor of about 20, which results in a significant reduction in fine-tuning. We have checked that the 
fine-tuning can indeed be mild by computing the logarithmic derivatives of the EW scale w.r.t. the 
microscopic parameters, in the framework of the 1-loop effective potential]^ Remember also that, as 
mentioned earlier, the Higgs potential is much less sensitive to the (Dirac) gluino mass in this model, 
compared to the (Majorana) gluino mass in the MSSM. 

It is also interesting to note that a combined global fit to Higgs properties, Br(i? — )■ ^^7), and 
the VT-mass, show a preference for 400-500 GeV degenerate stops, provided there is an additional 
mechanism to obtain a ~ 125 GeV Higgs eigenstate ^8|. This is precisely the situation in the bench- 
mark examples in Table [6} where the triplet provides the additional contribution to the Higgs mass 
(its heaviness being motivated by EWPT). Finally, note that since the couplings {Af , A^ grow with 
energy and since A^ is close to unity in Table |6} one find^ a Landau pole at around 10^ GeV (10^ 
GeV) for Benchmark I (Benchmark II), implying that additional new physics has to come in around 
those scales. This is consistent with our approach specified in the Introduction; we are primarily 
interested in understanding the nature of physics affecting the LHC, and are agnostic about effects at 
higher energy scales. Presumably a microscopic understanding of supersymmetry breaking within this 



^^We have used the fuh Coleman- Weinberg one-loop effective potential to compute the lightest CP-even Higgs mass 
eigenvalue. 
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^^We will present our results in more detail in 

^^We have computed the RGB's by implementing the full model in SARAH 59 60 
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setupp^ will provide insights into the nature of physics at such scales. 

Before moving on to discussing various aspects of collider phenomenology, it is worth commenting 
on the properties of the ~125 GeV eigenstate within the model. For the two benchmarks, it can be 
readily checked that this state is primarily a combination of and z/(a) with only a negligible Sr 
component. Thus, this state has properties very similar to the lightest CP-even state arising within 
the MSSM (with replaced by of course). As mentioned at the beginning of the section, the 
couphngs, production, and decays of this and the other scalar electroweak states will be studied in 



detail in another work 54 



7 Phenomenology 

In this section, we discuss several phenomenological features of the framework in which lepton number 
is related to the i?-symmetry {R = Ri). Since the phenomenology of this class of models is rather 
novel, in this paper we only outline the broad phenomenological consequences of the framework for 
collider and neutrino physics. A detailed treatment of these issues is provided in a companion paper |T], 
which studies the existing collider constraints on this class of models, as well as the various interesting 
signals which could be probed in the near future. 

The main qualitative features of the phenomenology of this class of models which sets it apart from 
traditional supersymmetric models like the (i?-parity conserving) MSSM (or many of its cousins like 



the NMSSM or models with extra vector-like mattei ), are the following: 



The existence of a "Dirac" structure in the gauge sector of the model. Among other things, 
this gives rise to a suppression of the production of squark pairs at the LHC compared to the 
Majorana case, which helps in relaxing the bounds on superpartners from current searches [l]. 
Another straightforward consequence of this is a suppression of those signals which depend on 
the Majorana nature of gluinos, e.g. same-sign (SS) dileptons. 

The i?-symmetry dictates a specific set of operators in the superpotential and Aoft, distinct from 
the "standard" cases. In particular, in the R = Ri realization, there exist "RPV" operators of the 
type XijkLiLjEl and X'^j^^LiQjDl in the superpotential, and -B/i^^ H^Li in the soft Lagrangian. 
Since these operators are consistent with the i?-symmetry (hence with lepton number), they 
cannot generate neutrino masses. Therefore, the sneutrino can have a significant vev in these 
models, thereby acting as a genuine Higgs field, in stark contrast to standard RPV models. 

Also, since the usual trilinear terms involving squark, slepton and Higgs fields are forbidden, 
there is no left-right mixing in the squark and slepton scalar mass-squared matrices. 

The existence of a sizeable (electron) sneutrino vev implies that some of the A and A' couplings 
are the lepton and down-type Yukawa couplings, respectively (which are well known up to tan /3). 
As explained in Sections [3] and |5} this implies that the flavor structure of the A and A' couplings. 



2.1 



At present, wc have only done a spurion analysis of supersymmetry breaking, in Section 
These are some of the popular models which could also give rise to Higgs mass near 125 GeV without much tuning. 
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as well as the various indirect constraints on these, are rather specific compared to standard 
RPV models. 

Furthermore, it implies that there is mixing between the neutrino (s) and neutralinos, and between 
the electron and the charginos, but with a different dependence on the parameters compared to 
that in standard RPV models, as explained in Section [3j 

A rich and interesting pattern of signatures results from such a structure. For example, decays of 
the "LSP"ghke X°+ —7- Zz/g, /ii^e; W e^, and — ?■ Tj^ T^e^iibR are prompt and have a significant 
branching ratio, unlike standard RPV models. 



The existence of an i^-symmetrj implies the conservation of two charges, the electric charge 



and the i?-charge, even after electroweak symmetry breaking. In particular, the neutralinos and 



charginos are Dirac in nature, and their interactions must conserve both charges (see Sections 4.1 



and 4.2). This gives rise to a rich and interesting pattern of decays of the neutralinos and 



charginos, and sleptons and squarks, which is different from that in the MSSM. 

In principle, flavor physics can be quite rich as well since bounds from flavor-violating processes 



are quite relaxed with an i?-symmetry 29 . However, we will not consider this in detail in this 
work. 

Although different subsets of the above set of signals can be mimicked by other models, the entire set 
of signals is rather unique. Hence, if the model is correct, it should be possible to distinguish this class 



of models from other models in the near future (more about this in Sections 7.3 and 7.5). 

Before going into more details about the phenomenology, it is useful to have an understanding of 
the spectrum of the model. Since the motivation is to build an electroweak-natural model, we consider 
a situation in which the third generation squarks are light (< 500 GeV). In fact, as will be shown 
in [l], the bounds on third generation squarks are weaker than 500 GeV, while that on the flrst and 
second generation squarks also turn out to be mild - in the range 600-700 GeV. This bound assumes 
a heavy Dirac gluino {M^ ^ 2 TeV), which is still consistent with naturalness due to the existence 



of a "supersoft" structure in the gauge interactions of the model 28 . The Dirac wino is expected to 
be heavy (> TeV) to satisfy electroweak precision measurements of the coupling of the Z to charged 
leptons for a reasonably large range of the sneutrino vev (see Section |5]). For concreteness we take 
M2 — 1.5 TeV. There are no direct bounds on the Dirac bino; however anticipating that the origin of 
its mass is tied to those of the wino and gluino, we take ~ 1 TeV for concreteness Since the /i 



parameter is directly connected to naturalness, we take it around the EW scale, jj, ~ 200 GeV. This 



'^^ "LSP" here stands for the lightest non SM-like superpartner which is charged under the SM. The "LSP" is really 
unstable, just as in RPV models. Note that the qualification of being charged under the SM is relevant because the 
gravitino can be the lightest BSM particle in many cases; however, in our framework, final states that include a gravitino 



have a negligible branching fraction and play no role in collider physics, unlike that in gauge mediation (see Section 7.5 1. 
Hence, we will reserve the term "LSP" for the lightest non SM-like superpartner charged under the SM, such as a 
neutralino or stau. This is phenomenologically useful since the "LSP" is the last step of the SUSY decay chains before 
producing a pure SM final state. 

^^It is assumed to be (explicitly) broken only by a very small amount, so for collider purposes the symmetry is exact. 

^^In Sectionlerj we have taken Mf = 200 GeV, but a CP-even Higgs near 125 GeV is also possible with A/f = 1 TeV. 
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implies that the hghtest non SM-hk^^ charginos and neutrahnos are mostly Higgsino-like. Finally, 
the sleptons are expected to be among the lightest particles in the BSM spectrum because of the 
close connection of the slepton sector with EWSB in our framework. Furthermore, a good degree of 
degeneracy among the three generations of sleptons is expected. Since the electron sneutrino plays 
the role of the down-type Higgs, electroweak naturalness requires its soft mass to be close to the 
electroweak scale. For concreteness, we take m| ~ m'i ~ (200 — 300 GeV)^ for all three generations. 
Thus, the lightest BSM particles consist of the sleptons, sneutrinos and Higgsino-like neutralinos and 
charginos. Depending on the situation, either the sleptons/sneutrinos or the Higgsino-like neutralino 
could be the "LSP". 



7.1 Summary of Bounds from Existing Searches 

In Ref. [l] we perform a more detailed analysis of the implications for the LHC of the leptonic f/(l)ij 
symmetry. We highlight here a subset of those findings. In the LHC context, the small i?-violating 
effects can be neglected, and the physics effectively exhibits a new conserved quantum number, the R- 
charge. As described above, we focus on the region of parameter space where the Dirac gaugino masses 
are in the 1 — 2 TeV range, which effectively decouples (in a first approximation) the bino-singlino, 
wino-tripletino and gluino-octetino Dirac states (though the latter can have some impact on squark 
pair-production through certain t-channel diagrams). As a result the lightest (higgsino-like) neutralino 
and chargino states play an important role, through the decays listed in the third item above. Such 
decay channels are also intimately related to the fact that down-type fermions get their masses from 
a sneutrino vev. We note that, due to its higgsino-like nature and the typical branching fractions, the 



LEP bounds on charginos are still stronger than those available at the LHC 61 , in our framework. 

Let us start by summarizing the results of the interpretation of the current ATLAS and CMS 
searches on the first and second generation squarks within the leptonic f/(l)_R model. The first point 



to notice is that, as emphasized in 22,23], the Dirac nature of the gluino results in a suppression of 



the squark pair-production cross section,^ By allowing for lighter squarks than in the MSSM, there is 
an important secondary consequence that compounds this effect. Namely, that the efficiencies of the 
current searches, which are optimized for MSSM-like production cross sections, can be significantly 
reduced due to the aggressiveness of the present cuts. The upshot is that squarks as light as 600 — 
700 GeV can be consistent with the various generic SUSY searches (involving jets, varying number of 
leptons and /E^), to be compared to the current limit of 1.4 TeV in the MSSM (when the gluinos are 
heavy) [62j . 

We also find that within our framework the bulk of the processes end up producing some amount of 
missing energy (in the form of neutrinos) , so that most of the present search strategies can apply with 
minor modifications. However, there are some topologies involving only visible particles (for instance 

^^Note that since the neutrino(s) and the electron mix with the neutrahnos and charginos respectively, technically 
these are the lightest neutralino(s) and chargino, respectively. 

^^The Dirac gluino pair-production cross section is larger than in the MSSM, but for decoupled squarks the current 
bounds on the gluino mass are only slightly stronger than for the Majorana case due to the steeply falling cross section 
with the Dirac gluino mass. In any case, we are taking gluinos at around 2 TeV to emphasize that a spectrum with 
squarks lighter than gluinos is fairly natural in our framework. 
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when the neutrahnos decay through their We channel, with a hadronically decaying W). It would 
be interesting to design search strategies for such "no ^t" channels. Another important example of 
no missing energy channel is the lepto-quark one. Interestingly, the RH strange squark has a sizable 
branching fraction into sr — )■ e^j. Current LQ searches have only a slightly smaller reach than the 
generic SUSY searches above (the latter assume degenerate squark masses), and a discovery in such a 
channel represents an exciting possibility. Besides allowing for a measurement of the LQ mass, such 
decay modes may allow to extract additional information from "mixed topologies", where one of the 
pair produced particles decays visibly through their lepto-quark mode, while the other one decays 
into states involving missing energy, which can be used for triggering and background reduction. One 



example of this sort is described in Subsection |7.4[ and illustrates how one may infer that the sneutrino 
vev is indeed large, a characteristic of the leptonic /^-symmetry. Other examples are presented in [l]. 

The third generation is likely of central importance for understanding the physics of EWSB. As 
was emphasized in Section [6T , the stops can be reasonably expected to lie in the few hundred GeV 
range, based on naturalness considerations, since the bulk of the lightest CP-even Higgs mass arises 
instead from the radiative corrections of a heavy triplet (with moderate associated fine-tuning). As 
it turns out, the lepto-quark nature of t^, bi and offers a powerful handle into this sector. As 
further explained in the companion paper [l] , depending on the sneutrino vev, such searches can easily 
cover the interesting expected range. Still, at present, the strongest constraint arises from the CMS 
direct sbottom search 63 , which can be interpreted as leading to a sneutrino- f ef -dependent bound 
on bi that varies between ~ 300 and 500 GeV, for representative values of the model parameters. 
The same search results in a bound on the RH sbottom of close to 500 GeV. The LH sbottom bound 
implies indirectly a lower bound on the LH stop a few tens of GeV larger, since due to the absence 
of LR mixing, the stop is always heavier than the sbottom. We therefore see that a large part of the 
interesting mass range has been tested, and that naturalness would lead us to conclude that a discovery 
should be possible in the relatively near future. We discuss a further interesting feature of a possible 



LQ signal in Section 7.3 One should also mention that most of the current dedicated searches for 
third generation squarks do not apply in a straightforward manner, but it should be possible to adapt 
them to cover more general possibilities than found in simple limits of the MSSM phenomenology. 



7.2 Resonant Slepton/Sneutrino Production 

One of the characteristic features of the presence of the X'^jf^LiQjDl operator is the resonant production 
of sleptons and sneutrinos, similar to that in RPV models. In this subsection, we would like to study 
the prospects for resonant slepton/sneutrino production at the LHC, assuming that the bounds on the 
A and A' couplings in Tables [3} |4] and |5] are saturated. As explained in Section |5| the sole exception is 
the X'^ii coupling, which we assume to be negligible (~ 0). Furthermore, for simplicity we will study 
the simple case where the sleptons and sneutrinos are the "LSP", so that the only decay modes of the 
slepton/sneutrino are via the A and A' couplings (resonant production occurs via A' couplings). Finally, 
we organize our analysis in terms of the two-body final states coming from slepton/sneutrino decays. 
For a given final state, we assume that only those A and A' couplings which could give rise to that 
particular final state are non-zero and that they saturate the bounds in Section [5j For example, for the 
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Figure 2: Upper bounds on <Jprod x BR for the e*/i^ (left panel) and 
electron sneutrino vev for two different values of the sneutrino masses 
i>T are produced, and the non- vanishing couplings are taken to be (A211, 



(right panel) final states as a function of the 



200 GeV and 300 GeV. For e^A*"^, 



both and 



A3]^]^, A322) for production and (A212, A231) 



for decay. For /i"*"/!", only is produced, and the non-vanishing couplings are taken to be A322) for production 

and (A322) for decay. The A and A' couplings identified with lepton and down-type Yukawa couplings are always present. 
The LHC bounds are shown as dashed lines. These bounds are provided by 64 for 



and by 65 for 



e /i^ final state, this would imply that the production of both z/^ and Ur is considered turning on those 
couplings which allow them to be produced (Ag^i, -^311 ' -^322)5 those which allow them to decay 
to e^/i^ (A2125A231). Of course, the A and A' couplings which are identified with the lepton Yukawa 
(Ai22,Ai33) and down-type Yukawa (A'^, A'122) -^133) couplings are always assumed to be present and 
non- vanishing ^ 

Comparing with the existing experimental bounds from the various two-body final states from the 
Tevatron and the LHC, we find that only the e^^"^ and fi~^fi~ final states provide constraints on the 
parameter space. Fig. [2] shows the bounds on the aprod x BR as a function of the (electron) sneutrino 
vev (varied within the allowed range) for the e^/x^ and fi^^~ final states. One finds that only values 
of the sneutrino vev close to the minimum value (maximum for tan /3) are allowed by the current LHC 
constraints. However, since these bounds are for values of couplings saturating the bounds in Section [5} 
this simply suggests that there are good detection prospects for these final states in the future, within 
this framework. 



7.3 Lepto-quark (LQ) Signals — i?-symmetry at the TeV Scale 

Another important consequence of the X^j/^LiQjDl operator is the presence of lepto-quark (LQ) signals. 
Again, since the LQD'^ operator is also present in standard RPV models, these signals are present in 
principle in RPV models as well. However, we will argue below that observation of certain LQ signals 
will in fact suggest the existence of an /^-symmetry in the TeV scale Lagrangian, which is, furthermore, 
tied to lepton number [R = Ri). This is made possible by a connection to neutrino physics, as will 

^^Note that we assume that if more than one sneutrino or slepton species can lead to the given final state, they have 
comparable masses but are split by an amount larger than their width, so that interference effects are negligible. This 
is well justified since the widths of the sneutrinos/sleptons are expected to be extremely small. 
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be explained. Finally, we will mention a number of situations in which further information from 
other channels would be required in order to rule out alternative interpretations. Fortunately, such 
information should be accessible at the LHC if the leptonic U{1)r symmetry is indeed at play. 

Because squarks have i?-charge 1 in the model, they also carry lepton-number since it is identified 
with the i?-symmetry {R = Ri). These scalar "lepto-quarks" are pair produced by QCD interactions, 
but can decay via the A' coupling above, thereby displaying their LQ nature. Such channels can be 
very important for the third generation^'' for two reasons - i) the bounds on A333 are quite weak from 
Section [sj and ii) the third generation Yukawa couplings are the largest, while the third generation 
squarks can naturally be the lightest. In particular, this means that and A333 can be and 
could give rise to the decays 

h bRU , Ir hiv , z/ = z/e, z/^ , (29) 

plus their conjugate processes. Since for an electroweak-natural model we expect the third generation 
squarks to have masses < 500 GeV, there are good prospects of observing these signals. The obser- 



vation of the LQ signal (29) may provide support for the existence of the U{1)r symmetry in the full 
Lagrangian at the TeV scale, which is a stronger conclusion than the one one could reach based on the 
observation of Dirac gluino signatures alone. The argument is the following. 



Let us start by assuming that some of the LQ signals (29) have been observed [either from Br or 



from the {tL,bL) doublet]. Although it may be hard to distinguish the decays in the second line in 



(29) from the "standard" SUSY decays involving a (massive) neutralino instead of the neutrino, the 
simultaneous observation of the (fully visible) decays involving a charged lepton may be taken, based 
on SU{2)l invariance, as indication that (at least part of) the missing energy signal is associated with 
a neutrino. We will also assume that the observed LQ is indeed a third generation squark (as opposed 
to first/second generation), and will further comment on this assumption below. 

Given the above LQ observation, it is possible to obtain ^2 

LR, 

a surprisingly large amount of information regarding the ^ ^ X ~ ^ 

structure of the soft supersymmetry breaking terms. This 6^ bji 

is because the LQD^ operators lead to a contribution to the nib ^ 

neutrino masses given by: ^ ^ ^ ^ ^ ^ ^ 

b"" b 

Am ^ 

^33)^ ^LR 7 = 19'? i"?n'i Figure 3: Contribution to the neutrino mass 

^ 167r2 m? ' - ' ' ' I J j^^j^ ^ gygY lepto-quark. 

where rnf^j^ is the left-right mixing in the sbottom mass-squared matrix. Note that in the /^-symmetric 
limit this left-right mixing is forbidden, but it is present in RPV models in general. For an electroweak- 
natural model, barring fine-tuned cancellations amongst several contributions, there are two ways in 



which Eq. (30) can be consistent with the upper bound on neutrino masses: i) a suppressed coupling 



^''^The RH strange squark can also display interesting LQ decay channels, see Jlj. 

^^Recall that when R = Ri, the couplings Xijj, j — 2,3, and X[i^f., k — 1,2,3, are identified with the lepton and 
down- type Yukawa couplings, respectively. 



26 



A^33) or a) a suppressed m1^. 

In RPV models, where typically = nib^Ah — /itan/3) ~ 0{mbmi), this would imply situation 
(i) - a very suppressed coupling A^33 (< 10~^). Here fi denotes the traditional "/i-term"|^ However, 
the observation of a LQ signal with a very suppressed A^33 coupling is only possible for a rather special 
situation, i.e. when the LQ decay channel (controlled by A'j33) of the relevant squark (sbottom here) 
has a significant branching ratio (BR) in spite of A^33 being very suppressed. We will discuss cases in 



Section 7.3.1 when this special situation is satisfied. In a typical situation, however, one expects at 
least some neutralinos and/or charginos to be lighter than the third generation squarks (which could 
have masses ~ 400 — 500 GeV). This is also motivated by electroweak-naturalness, where a natural 



solution of the EWSB minimization conditions typically requires /i < 200 — 300 GeV,^° implying the 
existence of at least one neutralino and one chargino lighter than the squark (sbottom here). Then 
standard SUSY two-body decays of the squarks are also open, and can be used to set a lower bound on 
A', based on the fact that the BR associated with the lepto-quark channel cannot be very suppressed 
if such a signal is seen at the LHC 



In Appendix [C| we estimate this lower bound on A'j33 for the two final states in (29 ) without missing 
energy - (i) charge |, fe_R, and (ii) charge -|, l^tR- For (i), we find A'j33 > 0.01 — 0.1, while for (ii) 

3' 'L 



the lower bounds are slightly smaller. Thus, for the charge |, Itbji final state, we estimate: 



A -/i tan/3 < (0.002 - 0.02) (^^^-^j GeV, (31) 

while the charge -|, l^tn final state could be consistent with a sbottom LR mixing term about two 
orders of magnitude larger. This implies that barring fine-tuned cancellations, both Ab and fi must be 



highly suppressed relative to ~ Msusy- Furthermore, note that the upper bounds (31) are valid 
at around the electroweak scale. Now, Ah in particular gets contributions from RG running from all 
Majorana gaugino masses as well as other A terms according to: 

dAh 1 f 32 „ , 14 , , , „ , 1 , , 

-gjM, + GgjM^ + -gfM, + 2y^At + Uy^A^ } , (32) 



dt 167r2 1 3 ^'^ " 15- 

where y and A (with appropriate subscripts) are the Yukawa couplings and A-terms,^^ respectively, 
while Mi are the SU{3)c x SU{2)l x f/(l)y Majorana gaugino masses. Even for a low scale cutoff 
A ~ 10 TeV, the individual contributions to Af, would be: 

AAb ~ {0.2 M3, 0.05 M2, 0.002 GeV , (from Majorana masses) , 



29 



More precisely, the relevant /i-term corresponds to the superpotential bilinear linking the Higgs doublets giving rise 



to up-type and down-type fermion masses, and is not the same as fi in (11). In fact, in i?-symmetric models, /t = 0. 
is the coefficient of the term HuRd in the superpotential in ( 11 ). 
^^The trilinear soft terms are defined with one power of the corresponding Yukawa coupling factored out, as is 
customary in the flavor diagonal case. 
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/ tan f3 



AAb ~ {0.04 A, 6 X 10"Mfe,4 x lO^M^} ( ) GeV , (from A terms) . (33) 



2 



The weaker bound in (31) would then imply 



/ 777," \ ^ 

(Ms.M^.Mi) < (0.1.0.4,10) (g^j-^j GeV. 
(A, A, A.) < (0.5.30.500) (^)' GeV . (34) 



which would become stronger by one order of magnitude if the stronger bound in (31) apphes. Thus 



in addition to upper bounds on Ab and p, from (31), we also get significant upper bounds on the three 



Majorana gaugino masses as well as on At. For larger tan /3 the bounds are even tighter and one 
could also infer an upper bound on Aj. as low as a few GeV. Since the absence of Majorana gaugino 
masses, A-terms, and the ft term are the hallmark of an i?-symmetry, the observation of a lepto-quark 
signal, via a connection to neutrino masses, allows one to build a strong case for the presence of an 
approximate t/(l)i? symmetry at the TeV scale. 



7.3.1 Special Cases 

As mentioned earlier, there exists special cases in which RPV models can also give rise to visible LQ 
signals, consistent with neutrino masses and other constraints. First, it may be possible for the LQ to 
correspond to a second generation squark, e.g. c^, decaying via the X'123 coupling as cl b^e']^ (and 
sl ^ bni^e)- The point is that such couplings are not constrained by neutrino mass bounds, and the 
constraints arise from very loose "single coupling bounds", or from "product coupling bounds" that 
can allow be sizeable if the other coupling is sufficiently suppressed. If couplings like X'123 are 

indeed sizeable, special kinematic configurations would not be required for the LQ channel to have a 
visible branching fraction. Nevertheless, observation of a LQ signal is only feasible if the LQ is well 
below 1 TeV. Thus, given the strong bounds on first two generation squarks in RPV models without 
an i?-symmetry, such an interpretation would require a significant mass splitting between cl and the 
other squarks. Such a situation would be distinguishable from the case where the LQ signal arises 
from third generation squarks. For instance, the "standard" decays of the LQ to the "LSP" (which 
decays further), for example cl ex? versus — )• txi, could distinguish between the two situations. 

Thus, we focus on the more natural SUSY interpretation of a LQ signal as a third generation 
squark, which might still be consistent with a very suppressed A'j33. This happens when all other decay 
modes of the squarks are suppressed, so that the LQ decay modes have a significant branching fraction 
even with suppressed A^gg couplings. Some situations in which this can arise are the following: 

• There is only a single LQ signal and it is the LSP. Then only LQ decay modes are available, if the 
LQ is bji. If the LQ is the doublet, one could hope to use the decay of the heavier into the lighter 
SU{2)l component, via an (off-shell) W, to extract additional information. However, besides the 
3-body phase space suppression factor, the corresponding partial width scales like (Am)^, where 
Am is the sphtting between the two SU{2) lepto-quark components. As a result, the derived 
lower bound on (A')^ is not useful to conclude that the LR mixing should be suppressed. 
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The observed LQ is 6/? while the doublet LQ's are too heavy, and/or the BR in their LQ channels 
is too small (so that they are not seen in those channels), and if in addition the LSP neutralino 
is almost pure wino or almost pure if^, then the 2-body decay of into the neutralino LSP is 
highly suppressed. This case could then be similar to the previous one, unless the second-lightest 
neutralino is lighter than b^. 

All neutralinos and charginos are heavier than the LQ (or very near threshold), even if the LQ 
is not the LSP. One then has to compare against loop-induced 2-body decays, or 3-body decays. 
Such decays are sufficiently suppressed that a sizeable BR in the LQ channel is allowed, consistent 
with neutrino bounds and without LR suppression. 

Even though the LSP is a neutralino, the LQ happens to be a (highly mixed) stop with the tx^ 
channel kinematically closed. One is then again left with loop-induced 2-body decays, or 3-body 



decays that can compete against the LQ signal even for rather small A' 66 . 

In order to further discriminate between the f/(l)ij-symmetry interpretation and the above possibilities, 
further information regarding the SUSY spectrum would be required. For example, the observation of 
the prompt decay of the neutralino to say We, which indicates the presence of an appreciable sneutrino 
vev, could be used to reconstruct the mass of the neutralino as well as shed light into the structure 
of neutralino interactions, providing support for this class of models over traditional RPV models. A 
thorough discussion of these issues is left for future work. 



7.4 LHC Signals of a Large Sneutrino vev 

We have seen in previous subsections that the model produces several distinctive signals at the LHC. 
Although standard RPV models can give rise to many of these signals in principle, we saw in Sec- 



tions 7.2 and 7.3 that the full pattern of signals is generically different. Furthermore, there is a 
particular signal topology which clearly distinguishes this model from standard RPV models. Not sur- 
prisingly, this difference arises due to the presence of a significant sneutrino vev, which is a distinctive 
feature of the model. 

Such a sneutrino vev opens up the neutralino decay modes into We, Zv and hv. As will be explained 
in detail in the companion paper jl], these decay modes have sizeable branching fractions when 
is the "LSP". We therefore focus on this particular situation. The decay mode into We points 
unambiguously toward a mixing between the electron and the charginos. Although this signal could 
be interpreted also in the context of a standard RPV model with the left-handed sneutrino acquiring 
a vev, since the bounds from neutrino masses on the sneutrino vev are so stringent in the standard 
case, the neutralino typically decays outside the detector or through a displaced vertex. Therefore, a 
prompt decay of the neutralino into We is a clear sign of a sizable sneutrino vev and therefore a hint 
of a leptonic i?-symmetry {R = Ri). 

Furthermore, even in the case where the decay X^^ — )■ Zu^ or hue is dominant (with the We 
channel suppressed, as can happen in some regions of parameter space), the observation of a "mixed 



topology" signa from the pair-production of third generation squarks, could provide a large amount 
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That is, a LQ decay channel on one cascade leg and a "standard" SUSY decay on the other one. 
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of information. For example if one of the pair-produced t^s decays via a LQ channe as in Section 7.3 
while the other one decays via tX°^ — t- t{Zi'e, hue}, this could be argued to provide evidence 

for a leptonic /^-symmetry. This is because an observation of such a signal would allow us to draw 
several conclusions: 

• First, there is a neutralino lighter than t^. In addition, the A' coupling of the LQ decay channel 
is large enough to give an observable signal. In particular, one can conclude that the magnitude 
of A' is comparable to that of the electroweak gauge couplings. 

• The invisible particle is most probably a neutrino. It cannot be a neutralino LSP (in an RPV- 
MSSM scenario) since it would have decayed promptly through the A' coupling above into b- 
quarks and a neutrino (but we are imagining that the Z ot h above have been reconstructed) . The 
invisible particle cannot be the gravitino either since the three body decay mode of the neutralino 
into bbu via a A' comparable to electroweak gauge couplings will dominate over the two body decay 

into GZ or Gh. For example, for a bino-like NLSP, r(x? ^ GZ) = "1^,, (l - ^) sin^ Ow, 
while T{xi hhu) ^ \'''^g''^ 12288 ^^mj "^^ere mi is the mass of the (off-shell) sbottom squark. For 

h 

X' comparable to electroweak gauge couplings and less than a few TeV, the decay width into 
the gravitino is always suppressed compared to the three body decay. 

• Finally, as argued earlier, an observation of the prompt decays Xi ~^ ^'^^ hv is strongly disfavored 
in standard RPV models from bounds on neutrino masses. 

Therefore, an observation of the above mixed topology would also strongly point towards a large 
sneutrino vev, and hence towards a leptonic i?-symmetry {R = Ri). 

7.5 Distinguishing from Other Models 

We now briefly discuss how some other popular supersymmetric models can be distinguished from this 
model. The "classic" and well studied cases of the "constrained MSSM (cMSSM)" and "phenomeno- 
logical MSSM (pMSSM)" are easy to distinguish from this model. This is because the presence of 
i?-symmetry in both the gauge and matter sectors, implies many signals which are quite different from 
the cMSSM and pMSSM models, such as: 

(a) Different production rates and decay modes for squarks. 

(b) Dirac gauginos with suppression of SS dileptons signals. 

(c) Resonant slepton/sneutrino production. 

(d) Lepto-quark (LQ) signals. 

(e) Decay of the "LSP" giving rise to fewer channels with^T, etc. 

^•^This has a significant branching ratio, especiahy if A333 saturates the bounds in Section [sj 
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The issue of distinguishing GMSB models from this model is more interesting. In "standard" GMSB 
models (without any i?-symmetry), signals (a)-(d) above should still allow us to easily distinguish it 
from this model. In /^-symmetric gauge mediation models (but with the /^-symmetry not related to 
lepton number, i.e. R ^ Ri), the signals (a) & (b) are the same as in this model. Signals (c) & (d), 
however, arise from the LQD'^ operators which are only present when R = Ri; hence signals (c) and 
(d) can be used to distinguish among the two models. 

Signal (e) deserves some more comments. In GMSB, the gravitino is the LSP, hence the NLSP 
(such as a neutralino or a stau) can decay to the gravitino via x? ZG, hG and f — )■ tG 67 . This can 
resemble some of the decays of the "LSP" in our model: X^~^ Zue, hue and — )■ T^z^e, as explained 
at the beginning of Section [rj However, within our framework, the decay modes X'^'^ — )• W~e~^, 

— i- tLbfi, etc., can have a non-trivial branching fraction |jlj, and do not give any when the W's 
decay fully hadronically, thus allowing for a full reconstruction. Therefore, these decays can clearly 
distinguish our model from GMSB models. Note that the gravitino is, of course, also present within 
our model, and is also very light. However, due to its suppressed couplings with matter, the branching 
ratios of the "LSP" to final states with a gravitino are extremely small and can be neglected- 
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7.6 The Case "5 = i?" 

We briefly mention some phenomenological features of the case R = R2, i.e. in which the baryon 
number is identified with an /^-symmetry. As mentioned earlier, this case has already been studied 
in 



26 and some aspects of the collider phenomenology have been studied in detail in 68 



In these models, the baryonic RPV operators X'ljf^U^DjD'^ are consistent with the /^-symmetry. 
Therefore, the "LSP" decays into jets via the A" couplings in this case, giving rise to final states with 
very small missing energy. Thus, they may naturally evade most of the current LHC bounds, and 
have generated renewed interest. However, it is important to note that there is a subtlety. In order 
to evade the bounds, the decay of the "LSP" has to be prompt, which puts a lower bound on the 



A" couplings. A" > 10 ^ — 10 ^ 21 . However, the indirect bounds on many of these couplings from 



neutron-antineutron (n — n) oscillations are orders of magnitude stronger than this 21 '^^ Thus, having 



a large enough coupling requires either a specific flavor ansatz [69j, or some other mechanism which 



relaxes the bounds, for example if i?-parity is broken collectively 70 . 

Identifying the baryon number with an i?-symmetry provides an elegant way to evade these bounds 
(this point was already made in |26]). Indeed, analogous to the L = R case, here the U^D^D'^ operator 
is consistent with the i?-symmetry which is identified with the baryon number; hence there are no 
bounds on A" couplings from baryon number violating processes like n — n oscillations. The baryonic 
i?-symmetry will be broken by the gravitino mass m3/2, hence n — n oscillations will constrain m3/2 
depending on the details of /^-breaking and mediation. For example, the contribution to n—n oscillation 



^"'Note that this is in contrast to standard GMSB models, and even to i?-symnietric GMSB models with R = Rq, 
where the only available decay mode of the NLSP is to final states with a gravitino. 
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This is true if the couplings are "flavor-generic" . 
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from tree level sbottom/gluino exchange (see Fig. |4]) provides an upper bound on A'/^^g 21 : 
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where Mg is the Majorana gluino mass and 
rriLji is the left-right mixing among the sbot- 
tom squarks. Thus, an C(l) A" coupling is pos- 
sible for Mg,mLR ~ few GeV. Within generic 
gravity mediation, this implies ^3/2 ~ few GeV 
while for anomaly mediation 7713/2 ~ 100 GeV. If 
m3/2 > mproton, then there is also no constraint 
from proton decay (to the gravitino, since in that 
case it would be the lightest baryon!). In Fig. |4] 
we also show how the n ~ n oscillation process is 
related to the single resonant production of hR. 
Such a diquark signal would formally play the Figure 4: Upper panel: n - n oscillations from sbottom 
analogous role of the lepto-quark signals in the exchange via A"i3, showing the Majorana gluino mass and 
L = R reahzation discussed in previous sections. insertions (as boxes). The arrows indicate the R num- 




ber flow. Lower panel: single sbottom resonant production 
via A'/i3. 



It would be interesting to try to use this feature 
to infer that the TeV scale Lagrangian does in- 
deed display an approximate i?-symmetry, as was done in the L = R case. This is significantly harder 
here, since the n — n constraint is intimately connected to the first generation, while establishing that 
the jets in a diquark signal are in fact connected to the first rather than the second generation is 
not quite feasible. Of course, even observing such a resonance is significantly more challenging than 
observing a LQ signal, especially below the 1 TeV scale (most recent LHC di-jet searches impose cuts 
above 1 TeV to control the QCD backgrounds). In this sense, a diquark signal arising from the first 
two generation squarks (which can be heavier consistent with electroweak naturalness), through A'/j^^ 
or \ii2-, can be interesting. However, since the squarks could well lie below 1 TeV in this model, the 
low mass region should be kept in mind for such searches. This would also apply for third generation 
squarks decaying into dijets. 



8 Dark Matter 

Since in this class of models the "LSP" decays, there is no WIMP Dark Matter (DM) candidate, 
similar to generic RPV models. However, other mechanisms for DM generation are available. We 
leave a detailed exploration of these issues for future work, and only briefly mention a few possibilities. 

First, even in the minimal setup considered in this paper, the gravitino can provide a natural DM 
candidate, as follows. Although we have not studied effects of i?-breaking in this paper in detail, 
such operators are necessarily present, as explained in Appendix |Aj One of the crucial consequences 
of these i?-breaking operators is the generation of neutrino masses as mentioned earlier. Since the 
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scale of i?-breaking is ultimately tied to m3/2, the upper bound on neutrino masses places an upper 
bound on m3/2, the precise magnitude of which depends on the details of i?-breaking mediation (see 
Appendix |A] for a discussion on two natural possibilities - i) generic gravity mediation, and ii) anomaly 



mediation). Ref. 25 has studied fitting the entire pattern of experimentally measured neutrino masses 
and mixing angles within this framework, and finds that 777.3/2 ^ 0{keV) for generic gravity mediation, 
and 7773/2 ^ 0{1 — 100) MeV for anomaly mediation. Assuming that the initial gravitino abundance 
is negligible and is such that it never reaches thermal equilibrium]^ gravitino masses in the keV-MeV 
range can indeed provide the DM abundance of the Universe by one of the following two processes: 



Thermal scattering of superpartners in the early Universe, as originally explained in [71j . In this 
case, the gravitino abundance depends linearly on the reheat temperature of the Universe Tr. 

Decays of superpartners, which are still in thermal equilibrium, to gravitinos - also known as 
Freeze-in (FT) ff^. In this case, the gravitino abundance is independent of Tr. 



Depending on the superpartner spectrum, one or the other process may dominate or they may both 
be comparable. Gravitino FI has been studied in [73] (for other discussions of the FI mechanism, 
see 74-76]). An important point to remember is that since technically the neutrino(s) are the lightest 



neutralinos in our model (by virtue of lepton number being an i?-symmetry, see Section 4.2), the 
gravitino is not absolutely stable. The dominant decay mode of the gravitino G is the process: G — )■ 
7 + z/, and it can be shown that the gravitino lifetime is sufficiently long so as to satisfy all observational 
constraints^ [77 . 



Finally, it is worth commenting on the tentative 7-line signal at around 130 GeV from the Galactic 



Center (GC) observed by many groups 78 81] in the FERMI-LAT data. It is clear that if the signal 
turns out to be correct (confirmed by FERMI-LAT), and if it is attributed to DM annihilation, then 
this cannot be explained within the framework above, at least in its minimal incarnation. Explaining 
the signal from DM annihilation within this framework would presumably require the existence of an 
appropriate additional (dark) sector. Studying these issues is left for future work. 



9 Conclusions and Future Directions 

Without a doubt, we have entered a data-rich era that is expected to finally unravel the mystery of 
electroweak symmetry breaking, i.e. uncover the physical microscopic mechanism underlying this well- 
established phenomenon. Although the discovery of a Higgs-like particle near 125 GeV and the absence 
of new physics so far provides some support for electroweak-tuned theories, it is still rather premature 
to abandon electroweak-naturalness. Indeed, Nature could be cleverly realizing an electroweak-natural 
model which manifests itself at the LHC in non-standard ways. In this work, we have studied one 
such elegant model, the defining feature of which is the existence of a continuous i?-symmetry, that 
coincides with lepton number when restricted to the SM sector. An important consequence of this is 

^^These conditions can be easily satisfied since the gravitino has extremely suppressed couplings. 

^^We disagree with earlier results for the decay width of the process G — > 7 + z^. We believe this discrepancy arises 
due to the earlier works not taking into account gauge invariance properly, in particular due to effectively using a 
non-gauge-invariant regulator. 
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that one of the left-handed sneutrinos gets a significant vev, which is not constrained by the upper 
bound on neutrino masses (we take the sneutrino getting a vev to be of the electron type, so that 
the corresponding charged lepton can easily get a mass from suppressed operators, while various other 
constraints are also satisfied). As a result, a large region of parameter space of the model is still viable, 
and leads to a rather rich and interesting phenomenology. 

The most important features of such a framework include i) "Dirac" gauginos (especially gluinos) 
due to "/^-symmetry in the gauge sector" , ii) Absence of certain scalar trilinear terms and the "stan- 
dard" /i-parameter (that we have called /t in this work) due to the "i?-symmetry in the matter sector" , 
iii) Existence of RPV operators of the type XijkLiLjE'^, X'-j^^LiQjDl and Bfi^l^HuLi consistent with 
the leptonic /^-symmetry R = Ri. Due to a significant sneutrino vev, the sneutrino acts as a Higgs 
field providing the down-type fermion masses through subsets of these operators, iv) Mixing between 
neutralinos and neutrinos, and between charginos and charged leptons. These features combine to 
give a rich and complex pattern of signals at the LHC, which is studied in a companion paper [l]. 
Although subsets of these features are shared by other models, the entire pattern of signals is rather 
unique and should be distinguishable from other models. Here we have highlighted how the (perhaps 
imminent) observation of lepto-quark signals could be construed as a powerful indication that the 
TeV scale Lagrangian indeed displays an approximate U{1)r symmetry of the leptonic type. We have 
also discussed how certain topologies at the LHC could be used to infer that the sneutrino vev is 
non- vanishing (and large), and also that a missing energy observation should be interpreted as being 
associated with a neutrino (as opposed to a neutralino or a gravitino). Further exploration of these 
issues will certainly be welcome in the near future. 

We emphasize that these models can easily accommodate the observed ~ 125 GeV Higgs-like signal 
with significantly less tuning than in the MSSM. The main ingredient is the existence of a scalar triplet 
which must be somewhat heavy (> 1 TeV), as required by current constraints on the p-parameter. 
If such a triplet scalar has an order one coupling to the Higgs, it can contribute significantly to the 
Higgs mass at the loop level (the i?-symmetry forbids other tree-level contributions). This is similar 
in spirit to using the stops for such effect within the MSSM, but with significantly less tuning than in 
that well-known case. This observation is more general than the particular model we are considering, 
but it fits rather nicely within the U{1)r framework. 

Finally, from the theoretical point of view, it would be very interesting to find a simple dynamical 
mechanism of supersymmetry breaking and mediation which preserves an approximate U{1)r symme- 
try, followed by an /^-symmetry breaking and mediation mechanism, which generates the structure of 
operators as envisioned in this class of models. 
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A i?-breaking Operators 



As explained in Section '2f2_. there are two natural ways in which i?-breaking can be transmitted to the 
visible sector: i) generic gravity-mediation, and ii) anomaly-mediation. 

Starting with generic gravity-mediation, we write arbitrary Planck suppressed couplings that are 
restricted only by the SM gauge symmetries, but which do not respect the (anomalous) U{1)r sym- 
metry (for concreteness, here we consider the case R = Ri). Nevertheless, we assume that <C Mp, 
so that the U{1)r symmetry is approximate in the observable sector, being dominated by the physics 
at M*. Thus, we have ?7(l)ij-violating superpotential terms: 



Wj^ = fiHuH, + fi[HM + ImsS^ + ^MtT^ + ^MoO^ + yf^H.Q.D''^ + y^^H,L,E^ + ^ ysS\ (36) 



which can schematically arise from 



/ "^'^ {^"^" + + S^ + T^ + 0' + -^^{HaQ^D] + if.L.E^ + 8')^ . (37) 

We see that ii[ ~ Fx/Mp ~ m3/2. In this case, we also see that the contribution to the "standard" 
(but f/(l)K-violating) down-type and lepton Yukawa couplings is yfj ~ y^j ~ m^i2/Mp < 10~^^ for 
^^3/2 ^ MeV. A similar suppression can be expected for the U^D'^D'^ superpotential operator and 
other i?-violating trilinears. There are also f/(l)i?- violating soft-breaking terms: 



VJ°" = B^iRH^R, + A^H^Q.U^^ + A%H,Q,D^^ + Al^H^UE': + A%,L,QjDl + A^^,L,L,E^, , (38) 
which can arise from 

/ "^'^ ^ {^"^-^ + ^ [^"^^^^ + ■■■]}' (39) 

giving contributions to the 6-terms of order rn'^^^j ^'^d very suppressed contributions to the A-terms, 
of order {m^i2/Mp)m^j2- However, there are also terms of the form 



/ 



d'^e ^ [H^Q.U^ + ■■■] , (40) 



that give contributions to the A-terms of order Again, we consider only Planck suppressed 

couphngs between the hidden and observable fields, so that no JcfOXHuHa operator, which would 
induce a too large Bfip term, exists. Similarly, there can exist f/(l)ij-violating gaugino Majorana 
masses ^M^A^A", induced by 

I Jd'O^y^'^yVa + i^.c. (41) 

so that Ma ~ ■^3/2- Thus, the scale of f/(l)ii;-preserving terms is taken to be MgusY, which is assumed 
to be much larger than the scale of t/(l)ij-violating operators, set by m3/2. All possible dimensionful 
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terms allowed by the gauge symmetries are allowed, and are induced at least at order 7713/2- Note 
that singlets in theories of the present type have been argued to be safe from radiatively-induced 
tadpoles 46 . Finally, as mentioned in the main text, the experimental upper bound on neutrino 



masses puts an upper bound on 7713/2 in this case to be 0(keV) 25 



In the case of anomaly-mediation, we imagine a situation in which the "tree-level" transmission 



of i?-breaking by operators (37), (39), (40) and (41) is very suppressed, e.g. due to sequestering 



Then anomaly-mediation generates essentially the same operators (40) and (41), except that they 



are generated at one-loop rather than at tree-level. So, the coefficients above are suppressed by the 



one-loop factor 



1 

167r2 ' 



giving rise to an associated mass scale: 



j^anomaly ^anomaly 



?7l3/2 
1671^ 



(42) 



In this case, the upper bound on 7713/2 from the neutrino mass constraint is 0{1 — 100) MeV [25]. 



B A Flavor Ansatz 



In Section 5^2, we saw that there exist a number of constraints on the A and A' couplings from various 
flavor-violating processes. In Section [7], we also studied phenomenological consequences of an ansatz 
in which all the A and A' couplings are assumed to saturate the bounds in Tables [3| |4] and [s] (except 
which is taken to be negligible.). However, it is also possible to consider other simple ansatze 
about the flavor dependence of these couplings, and we discuss such an option here which allows a nice 
understanding of the relative magnitudes of the various couplings and also satisfles existing constraints. 
Also, it has the advantage of further reducing the number of independent A and A' couplings. 



We imagine that the model described by the superpotential in (11) has a flavor symmetry Gp = 
SU (3)q X SU (3)[/ X SU (3) d x SU (3) l x SU (3) e at some high scale which is broken by scalar components 
of "flavon" superfields. This implies that the up-type Yukawa couplings and the A and A' couplings 
(which include down-type and lepton Yukawa couplings) should be thought of as flavon (super) flelds 
transforming in a non-trivial representation oi Gp- Of course, the symmetry group Gi? is broken when 
their scalar components get vev''s and give rise to the numerical values of the Yukawa couplings. The 
representations in Table [T] allow all the relevant operators in the superpotential in (11). 



In particular, the above choice gives rise to the following structure for the coefficient of LiLjE^ 
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^iifc _ ^ijn yk ^ ^43) 

where e^^^ stands for the antisymmetric tensor with three indices, and Yj^ transforms in the (□, □) of 
SU{3)l X SU{3)e- Then, using the SU{3)l and SU{3)e rotation freedoms for Lj and E'^, respectively, 

^^In this section only, we use subscripts (superscripts) to denote transformation in the fundamental (anti-fundamental) 
of SU{3). 
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Table 7: Representations of the fields and Yukawa spurions under Gp- 
it is possible to choose a convenient form of Y,}^. In particular, by choosing: 
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1/2 



n ■ 

n 

n 



1 , 

2 , 
3, 



(44) 



one gets a simple structure of the Xijk couplings, in which the matrix Xijk is a diagonal matrix corre- 



(e) _ 



sponding to the usual leptonic Yukawa couplings appearing in PVyukawa in Eq- (|15|) with ?/2 
and 2/3 
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(e) 
1/2 



Furthermore, the matrices X2jk and As^fc have only one independent Xijk coupling: 
A231 = — A321 = yv, the rest either vanish or are related to the previous Yukawa couplings. 

The LQD'^ operators are less constrained by these arguments. However, it is convenient to make 
the following simple ansatz: 



(A 



i\ijk 



y^ {Y.y" 



(45) 



where i/i are defined in (44), and (Y^)-^'^ transforming in the (□,□) of SU{3)q x SU{3)d- Again, by 
choosing an appropriate SU{3)q x SU (3) d rotation, (YdY'' can be brought into the form diag(yi, ^2, ^3), 
which results in a very simple structure of the A^ couplings as well. In particular, the matrix X\ 



is diagonal and corresponds to the usual down- type quark Yukawa couplings in WVukawa in (15). The 



remaining A' couplings are completely fixed by these down-type Yukawa couplings and the numbers 



?/2, 1/3} in (44 ). Since y2 and y^ are the muon and tau Yukawa couplings, the only other independent 
input to fix all the remaining ones is yi = A231 = — A321. 

Since the above ansatz provides a concrete determination of all the couplings, it is straightforward 



(but tedious) to check constraints from all combinations of A and A' couplings in 21,51-53 . We find 



that all such constraints are readily satisfied, except that coming from IA231 X'^nl, which is violated by a 



small amount. This, however, just means that (45) is only approximate and receives some corrections. 



In any case, the main motivation for the ansatz (45) is to simplify the structure of the A' couplings 



^^This can be achieved by diagonalizing 1^ by bi-unitary transformations, followed by the exchange of the second 
and third rows, with an appropriate sign. 

^°Note the electron mass arises not from superpotential couplings, but from SUSY breaking operators; see Section [i] 
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and perhaps suggest a possible high-energy rationale for their structure. 

Finally, since the only new coupling within this approach is yi, it is important to know what 
one would generically expect for its size. In fact, the bound on the product |A']^33A233| = yb|A233| < 
1.1 X 10"^ |53] puts a bound on yi, since in our framework ?/& A233 — ^^^^ — (1.53 x 10^^ sec^ (3)/yi. 
Therefore, the bound implies that the largest non- Yukawa A' coupling (A333) is roughly given by: 

A333 ~ ^ ~ 0-011 ^ ' (46) 

which is quite small. However, this bound is specific to the previous ansatz and, experimentally, 
A333 can be significantly larger, as described in the main text (where we do not use the ansatz of this 
section). We still describe it here as it is simple and will also be considered in our LHC phenomenology 
study in 111. 



C Lower Bound on A'33 given the Observation of an LQ Signal 

In this appendix we argue that typically the size of the relevant A' couplings cannot be extremely 
suppressed, if a LQ signal was observed at the LHC. In order to estimate a lower bound on (A')^ such 
that the BR in the LQ channel is not very suppressed (and could be observed in the near future), we 
explore several possibilities, depending on whether the lepto-quark is t^, bi or Br. To be definite we 
estimate the "minimum" value of (A')^ by requiring that the partial decay width in the lepto-quark 
channel, given by Flq = [(A')^/167r] m-LQ equals one of the standard R-parity conserving 2-body decay 
widths 82 , r2-body, interpreted within the MSSM: |^ 



2 1 QtT 

(•^min) = — r2_body • (47) 

We show in Fig. [5] the result for 500 GeV lepto-quarks, when comparing to their partial decay widths 
into a neutralino plus quark (left column) and into a chargino plus quark (middle column). We 
have performed a scan over Mi, M2, G [0,600] GeV and tan/3 G [3,50], diagonalizing the MSSM 
neutralino and chargino mass matrices to find the spectrum and composition of the eigenstates for each 
parameter point. We compare to the dominant neutralino or chargino channel, and plot the (A'^^j^)^ 
defined above, as a function of Am = m-LQ — m^, where m-LQ = 500 GeV is the lepto-quark mass and 
is the appropriate neutralino or chargino mass. We also show curves corresponding to the limiting 
cases in which the neutralino is pure bino, pure wino or pure Higgsino (/i„), and also when the chargino 
is pure gaugino or pure Higgsino. For each scanned parameter point we have also estimated {X'^^^^Y 
based on the largest partial decay width of any of the neutralino and chargino channels (shown in the 
right column plots). We will use the latter as our estimate for {\'^[^)'^- 

We see from the plots in the right column that in the case that the lepto-quark is a SU{2)l 
doublet, (Aj^jjj)^ is typically above 0.01 (unless the decay is very close to threshold). We also note that 

^^We show the results within the MSSM structure to illustrate what an interpretation outside the U{1)b, symmetric 
framework would entail. A similar analysis within the L — R model results in qualitatively similar features regarding 
the expected sizes of the A' couplings. 
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Figure 5: Required value of (A')^ such that the partial width in the lepto-quark signal equals the standard 2-body 
decay width into neutralinos/charginos for a 500 GeV (upper row), (middle row) and bn (lower row), within the 
RPV-MSSM. The third column uses the largest of the neutralino/chargino channels to bound (A')^ for any given point 
in the random scan over Mi, M2, /i and tan/3. The smooth curves indicate the limit where the neutralino is pure bino, 
pure wino and pure higgsino, as indicated, and similarly for the chargino being pure wino or higgsino. 



the bulk of the cases has an even larger (A^j^)^ > 0.1. Our estimate given in (31) is obtained by using 
^2 



0.01 - 0.1 in (30), together with Am^ 



0.1 eV. We also see from the lower row in Fig. |5j 
that when the lepto-quark is fo/j, our estimate for (A'^j^j^)^ is weaker, although still larger than what 
the neutrino bound requires if there was no suppression from LR mixing in the sbottom sector. These 
two cases could be distinguished by the type of lepto-quark signal: hi for the doublet (we expect the 
mass splitting between ti and 6l to not exceed a few tens of GeV), and tl for the singlet. Both cases 
should be accompanied by a 6 + I^t signal, that would suggest that the missing energy comes from a 
neutrino. 
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